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Abstract 


The  purpose  of  this  thesis  was  to  determine  (both  experimentally  and  analytically)  the 
initiation  and  progression  of  failure,  stress-strain  response,  and  the  failure  loads  of 
Graphite/Polyetheretherketone  (Gr/PEEK)  laminates,  incorporating  an  eccentric  0.4  inch 
circular  discontinuity,  loaded  in  axial  tension  at  room  temperature/)  The  ply  lay-ups  of  these 
specimens  were  [0°16],  [90°  jg],  [±45°]4S,  and  [0°/±45o/90o]2s^For  each'of  theseply  -  o  y 

jps;  three  values  of  eccentricity  were  considered  (the  three  values  of  eccentricity  were 
determined  by  the  hole  location  within  each  specimen).  In  addition,  experimentation  was 
conducted  to  study  the  effects  of  boundary  conditions  on  the  failure  characteristics  of  the 
Gr/PEEK  laminates;  this  was  accomplished  through  the  use  of  a  special  mounting  fixture 

which  allowed  in-plane  rotation  of  the  specimens.  Finally,  experimentation  was 

V  9—'- 

conducted,  using  photoelasticity,  to  verify  the  gross  stress  states  of  the  Gr/PEEK 

laminate/predicted  by  the  analytical  study.  N 

C  '  ”  - - " 

Analytically,  a  nonlinear  material  finite  element  program  was  used  to  predict  the 

initiation  and  progression  of  failure,  stress-strain  response,  and  the  failure  loads  of  the 

Gr/PEEK  laminates.  In  addition,  the  effects  of  boundary  conditions  on  the  failure  modes 

of  the  Gr/PEEK  laminatesjwas  studied  analytically.  And  finally,  the  gross^tress  states 


of  the  Gr/PEEK  laminates  were  considered  in  the  analytical  portion  of  this  thesis,  f  'C  ft  ) 

The  experimental  and  analytical  results  were  then  compared.  The  initiation  and 
progression  of  failure,  stress-strain  response,  and  the  failure  loads  of  the  (Gr/PEEK) 
laminates  compared  quite  well  for  the  [0°16],  [90°16],  and  [0°/±45o/90°]2s  laminates.  For 
the  [±45°]4S  laminates,  the  analytical  predictions  for  the  failure  loads  underestimated  (from 
16.6%  to  40.5%)  the  experimentally  obtained  failure  loads.  Furthermore,  it  was  observed 
(both  experimentally  and  analytically)  that  the  failure  modes  and  failure  loads  for  the 
Gr/PEEK  laminates  were  not  appreciably  affected  by  the  boundary  conditions  caused 
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through  the  use  of  the  mounting  fixture  (in-plane  rotation  allowed).  Finally,  there  was 
good  agreement  between  the  experimental  and  analytical  predictions  of  the  "gross"  stress 
states  of  the  Gr/PEEK  laminates. 

Since  eccentricity  effects  on  the  failure  modes  and  failure  loads  of  the  Gr/PEEK 
laminates  were  of  interest  in  this  study,  the  analytical  study  attempted  to  predict  the 
behavior  of  the  Gr/PEEK  laminates  due  to  an  eccentric  circular  discontinuity.  The  eccentric 
circular  discontinuity  created  a  centroid  (area)  shift  in  the  Gr/PEEK  laminates.  This  shift  in 
centroid,  combined  with  the  axial  loading,  created  a  bending  moment  which  varied 
depending  on  the  eccentricity  of  the  circular  discontinuity.  The  analytical  study  was  able  to 
accurately  predict  the  failure  modes  and  failure  loads  resulting  from  eccentricity  effects. 
Although  the  final  failure  mode,  and  failure  load,  was  available  from  experimentation,  the 
failure  progression  resulting  from  eccentricity  effects  was  not  observed  during  exper- 
mentation.  All  experimentation  was  recorded  on  video  tape  for  post  failure  analysis. 
However,  failure  was  too  rapid  to  observe  failure  progression  (by  viewing  the  video  tape 
frame-by-frame). 


A  STUDY  OF  FAILURE  CHARACTERISTICS  IN  THERMOPLASTIC  COMPOSITE 
LAMINATES  DUE  TO  AN  ECCENTRIC  CIRCULAR  DISCONTINUITY 


I.  Introduction 

In  weight-sensitive  applications,  such  as  aircraft  design,  accurate  strength  and  failure 
characteristics  must  be  ascertained  before  composite  materials  can  be  used  to  their  full 
advantage.  Therefore,  future  design  advancements  rely  heavily  upon  our  ability  to 
accurately  predict  strength  and  failure  characteristics  of  composite  materials. 

Fiber  reinforced  resin  composites  typically  have  high  strength-to- weight  and  stiffness-to- 
weight  ratios.  For  this  reason,  they  have  become  important  in  the  aerospace  industry.  As 
more  demands  are  made  upon  new  generation  aircraft,  we  must  continue  to  develop  new 
composite  materials  with  desirable  strength  characteristics.  To  this  end,  Gr/PEEK  has  been 
developed.  This  polyether  based  thermoplastic  material  has  shown  several  advantages  over 
graphite  epoxy  (Gr/Ep).  Among  these  advantages  are:  ease  of  fabrication,  lower  weight 
(per  part),  higher  operating  temperatures,  increased  resistance  to  delamination,  and  higher 
fracture  toughness  [9]. 

Past  research  on  Gr/PEEK  laminates,  containing  circular  discontinuities,  has  shown  that 
the  nonlinear  relationships  between  stress  and  strain  must  be  considered  [10,13,17]. 
Therefore,  a  nonlinear  material  finite  element  program  [23,25]  was  employed  in  this  study. 
A.  Purpj?S&. 

The  purpose  of  this  thesis  was  to  determine  (both  experimentally  and  analytically)  the 
initiation  and  progression  of  failure,  stress-strain  response,  and  the  failure  loads  of 
Gr/PEEK  laminates  containing  an  eccentric  0.4  inch  circular  discontinuity.  These 
Gr/PEEK  laminates  were  loaded  in  axial  tension  at  room  temperature.  The  ply  lay-ups  of 
the  laminates  were  [0°ig],  [90°ig],  [±45°]4S,  and  [0°/±45o/90°]2s-  For  each  of  these  ply 


lay-ups,  three  values  of  eccentricity  were  considered  (the  three  values  of  eccentricity  were 
determined  by  the  hole  location  within  each  specimen). 

The  experimental  and  analytical  objectives  of  this  study  were  as  follows: 

(1)  Investigate  the  initiation  and  progression  of  failure,  and  the  failure  loads  of  the 
Gr/PEEK  laminates. 

(2)  Investigate  the  boundary  condition  effects  on  the  failure  characteristics  of  the 
Gr/PEEK  laminates. 

(3)  Determine  the  "gross"  stress  states  of  the  Gr/PEEK  laminates. 

(4)  Compare  the  experimental  and  analytical  results. 

B.  Background  and  Overview. 

Composite  materials  contain  two  or  more  materials  and  have  been  in  use  for  centuries. 

Prehistoric  civilizations  used  straw  in  mud  bricks;  plywood  existed  in 

early  Egyptian  construction;  and  medieval  knights  used  combinations 

of  various  metals  in  their  swords  [13]. 

For  the  purposes  of  this  thesis,  a  composite  material  is  a  material  consisting  of 
continuous  reinforcing  fibers  in  a  matrix  material.  The  fibers  and  matrix  are  considered,  on 
a  macroscopic  scale,  to  form  a  useful  material  which  acts  as  a  single  structural  element. 

One  advantage  of  composite  materials  over  "typical"  isotropic  materials  (aluminum,  steel, 
etc.)  is  that  the  composite  material  can  be  optimized  (tailored)  to  meet  a  specific  appli¬ 
cation.  Therefore,  the  fibers  of  a  composite  material  are  typically  oriented  so  that  the 
directional  dependence  of  strength  and  stiffness  of  the  composite  matches  the  loading 
environment. 

In  recent  years,  the  term  "advanced  composite  material"  has  been  used  to  describe 
composites  with  very  high  strength  and  stiffness  fibers.  These  advanced  composites 
typically  use  a  plastic  (epoxy)  or  metal  matrix.  The  Gr/PEEK  used  in  this  study  is 
considered  to  be  an  advanced  composite  material.  Gr/PEEK  consists  of  continuous  AS-4 
graphite  fibers  in  a  polyetherttherketone  (PEEK)  matrix. 


Since  composite  materials  can  be  tailored  to  match  the  loading  environment,  they  usually 
have  high  strength  and  low  weight.  Therefore,  as  previously  mentioned,  composite 
materials  have  found  widespread  use  in  the  fabrication  of  component  parts  for  the 
aerospace  industry.  Examples  of  some  of  these  components  include  the  wings  of  the  X-29 
aircraft,  and  the  horizontal  and  vertical  stabilizers  of  the  F- 15  and  F-16  aircraft. 

These  aircraft  components  are  usually  connected  to  other  components  through  adhesive 
bonds,  mechanical  fastners  (such  as  bolts),  or  a  combination  of  adhesive  bonds  and 
mechanical  fastners.  Since  mechanical  fastners  require  fastner  holes,  successful  designs 
rely  heavily  upon  the  ability  to  analytically  predict  the  stress  distribution  around  these  holes 
which  result  from  service  loads.  The  analytical  predictions  of  load  redistribution  around  a 
fastner  hole,  in  composite  materials,  is  not  a  simple  task. 

Although  there  are  many  techniques  for  predicting  the  ultimate  strength  in  composites, 
in  general,  two  techniques  are  commonly  used.  These  include  a  fracture  mechanics 
approach  [3,12,28]  and  a  finite  element  (damage  zone)  approach  [6,21,22,23,25]. 

Several  commonly  used  fracture  models  for  predicting  the  strength  of  composites  con¬ 
taining  circular  holes  are  reviewed  by  [3].  In  this  reference,  "serious  questions  are  raised 
as  to  the  applicability  of  classical  fracture  mechanics  to  composites  [3]."  Fracture  models 
typically  do  not  address  the  failure  modes  of  composites,  and  as  discussed  by  [3]: 

Semi-empirical  fracture  models  which  attempt  to  predict  the  notched 
strength  of  composites  do  not  address,  but  rather  by-pass,  the  micro- 
and  macro-failures  associated  with  the  crack  extension  process  [3]. 

Furthermore,  typical  fracture  models  require  correction  factors  which  are  often  arbitrary 

and  subject  to  question  [20] .  Therefore,  as  previously  mentioned,  the  analytical  portion  of 

this  study  made  use  of  finite  element  (damage  zone)  methods. 

For  several  reasons,  the  nonlinear  material  finite  element  program  which  was  used  in 
this  thesis,  was  ideally  suited  to  study  the  Gr/PEEK  laminates  considered  in  this  study. 
First  the  nonlinear  stress-strain  response  of  the  Gr/PEEK  laminates  was  modeled  by  piece- 


wise  cubic  spline  interpolation  functions.  The  use  of  piece-wise  cubic  spline  interpolation 
functions  to  represent  the  experimental  stress-strain  data  for  Gr/PEEK  laminates  [13], 
provided  "smoothness"  (continuity  of  slope)  over  the  entire  range  of  the  stress-strain 
curves.  Secondly,  the  finite  element  program  (unlike  fracture  methods)  properly  modeled 
the  damage  accumulation  process  of  the  Gr/PEEK  laminates.  This  damage  accumulation 
(failure  progression)  was  modeled  as  areas  of  failure  (element  failures)  within  the  finite 
element  model.  And  finally,  the  geometry  of  the  specimens  including  the  physical 
dimensions,  boundary  conditions,  and  ply  orientations  of  the  laminates  were  all  easily 
modeled  through  finite  element  techniques.  The  numerical  portion  of  this  study  (including 
the  development  of  the  finite  element  models)  is  discussed  in  detail  in  Section  III. 

As  mentioned  previously,  numerous  failure  theories  are  available  for  use  with  finite 
element  methods.  Over  30  failure  theories  exist  for  composites,  with  another  12  theories 
used  to  predict  post-failure  behavior  of  composites  [10].  Some  of  these  theories  may  be 
used  only  in  special  cases,  and  some  are  not  valid  for  use  with  composite  materials  [10]. 
This  study  made  use  of  a  total  strain  energy  failure  criterion  [21,22,23,25]  which  was 
incorporated  in  the  nonlinear  material  finite  element  program.  This  total  strain  energy 
criterion  is  discussed  in  detail  in  Section  II D. 

In  addition  to  the  analytical  study,  experimentation  was  also  conducted  to  verify  the 
analytical  predictions  of  initiation  and  progression  of  failure,  stress-strain  response,  failure 
loads,  and  the  "gross"  stress  states  of  the  Gr/PEEK  laminates  studied.  The  experimental 
portion  of  this  study  is  described  in  Section  IV. 

Section  VI  presents  the  results  from  the  analytical  and  experimental  studies,  and  the 
conclusions  from  this  thesis  are  summarized  in  Section  VI. 
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This  thesis  made  use  of  the  following  theories:  mechanics  of  composite  materials,  linear 
finite  element  theory,  nonlinear  constitutive  relationships,  and  a  strain  energy  failure  theory 
for  composite  laminates  [21,22,23,25], 

The  mechanics  of  composite  materials  can  encompass  both  micromechanics  and  macro- 
mechanics.  This  thesis  was  primarily  directed  toward  the  macromechanical  behavior  of 
laminates.  However,  for  completeness,  the  micromechanics  of  composites,  as  it  applies  to 
this  thesis,  will  be  discussed. 

Linear  finite  element  theory  will  be  discussed  briefly  as  a  "building  block"  to  a  dis¬ 
cussion  of  the  nonlinear  relationships  between  stress  and  strain.  Finally,  a  discussion  of 
the  failure  criteria  and  the  unloading  options  utilized  by  the  nonlinear  finite  element  program 
will  be  discussed. 


A.  Mechanics  of  G 


ite  Materials. 


(1)  Micromcchanical  Behavior  of  Composite  Materials. 

The  term  micromechanics  is  often  used  to  describe  the  behavior  of  composites  wherein 
the  interaction  of  the  constituent  materials  is  considered  on  a  microscopic  scale.  However, 
since  this  thesis  utilized  a  macromechanical  approach,  only  the  applicable  portion  of  the 
theory  of  micromechanics  will  be  discussed.  In  particular,  the  transfer  of  load  between 
fibers  is  of  primary  concern. 

An  aromatic  polymer  composite  (APQ  is  a  composite  which  is  constructed  with  a 
polyaromatic  resin  combined  with  continuous  fibers  [13].  The  Gr/PEEK  used  in  this  study 
was  APC-2  and  consisted  of  a  polyetheretherketone  (PEEK)  matrix  containing  Hercules 
AS-4  graphite  fibers.  For  Gr/PEEK  to  resist  externally  applied  loads,  the  fiber  and  matrix 
must  function  in  an  integrated  manner.  This  requires  that  the  load  be  transferred  from  fiber 


to  fiber  throughout  the  matrix  as  well  as  through  the  fiber-matrix  interface  [4,26].  In  a 
typical  composite  lamina  with  unidirectional  fibers,  the  fibers  are  surrounded  by  the 
adjacent  matrix  material  as  depicted  in  Figure  2-1.  The  fibers  in  Gr/PEEK  are  the  principal 
load-carrying  agent  The  matrix  serves  several  purposes  including  fiber  support,  protection 
of  fibers,  and  transmitting  load  between  fibers.  When  fiber  failure  occurs,  the  capability  of 
the  matrix  to  transfer  load  between  fibers  is  very  important  One  model  which  describes 
the  transfer  of  load  after  fiber  failure  is  depicted  in  Figure  2-2  [12].  In  this  model,  load 


Figure  2-1.  Lamina  with  Unidirectional  Fibers 


transfer  is  accomplished  by  shearing  stress  which  develops  in  the  matrix  material.  When 
composites  are.  subjected  to  high  temperatures  and/or  humidity,  this  interfacial  shear 
strength  can  be  reduced.  Therefore,  the  fiber-matrix  interface  is  often  times  considered  to 
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Figure  2-2.  Broken  Fiber  Effects 


be  the  weak  link  in  a  composite  lamina.  Numerous  theories  have  been  implemented  on 
how  load  transfer  between  fibers  is  achieved.  Two  of  these  theories  are  the  "shrink-fit" 
theory  and  "chemical  coupling"  theory,  both  of  which  are  discussed  by  [26]. 

Finally,  from  a  micromechanics  point  of  view,  the  properties  of  a  composite  depend 
upon  the  individual  properties  of  its  constituents.  Typically,  "the  composite  stress  and 
strain  are  defined  as  the  volume  averages  of  the  stress  and  strain  fields,  respectively,  within 


the  representative  volume  element  [27]."  Therefore,  it  is  necessary  to  find  relations 
between  the  composite  stresses  in  terms  of  die  properties  of  the  constituents.  A  major 
problem  with  micromechanics  is  the  inability  to  accurately  describe  the  boundary  conditions 
which  exist  between  each  constituent  within  the  material.  An  exact  analysis  based  on 
micromechanics  is  therefore  difficult  at  best.  For  this  reason,  this  thesis  made  use  of  the 
theory  of  macromechanics. 

(2)  Macromechanical  Behavior  of  Composite  Materials. 

Since  composite  materials  are  normally  nonhomogeneous  and  anisotropic,  the  theories 
and  formulas  with  which  we  have  dealt  with  isotropic  materials  cannot  generally  be 
utilized.  The  macromechanics  approach  concerns  itself  with  "gross”  composite  material 
behavior.  This  approach  does  not  recognize  the  individual  constituent  materials  and  the 
composite  material  is  considered  to  be  anisotropic. 

Since  most  practical  composites  are  uniformly  constructed 
(at  least  within  each  layer,  if  laminated),  each  layer  of  the 
material  can  usually  be  assumed  to  be  quasi-homogeneous 
on  the  macro  scale  [1], 

The  individual  layers  within  a  composite  laminate,  as  depicted  in  Figure  2-3,  are  bonded 
together  to  act  as  a  single  structural  element  Typically,  the  laminae  principal  material 
directions  arc  oriented  so  that  the  structure  is  capable  of  resisting  load  in  several  directions. 
The  stacking  sequence  of  laminae  can  therefore  be  used  to  "tailor  the  directional  dependence 
of  strength  and  stiffness  of  a  material  to  match  the  loading  environment  of  the  structural 
element  [12]."  The  complexity  arising  in  composite  laminates  is  the  introduction  of  shear 
stresses  between  layers.  These  interlaminar  shear  stresses  develop  due  to  the  tendency  of 
each  individual  ply  to  deform  according  to  its  respective  material  properties.  These 
shearing  stresses  are  largest  at  the  edges  of  the  laminate  and  often  lead  to  delamination. 
Additionally,  the  transverse  and  normal  stresses  which  aiise  from  uniaxial  loading  can  also 
lead  to  delamination  [12].  Therefore,  the  stacking  sequence  of  laminae  not  only  affects 
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Figure  2-3.  Laminate  Construction 


directional  dependence  of  strength,  but  also  affects  the  interlaminar  stresses  [28]. 

To  understand  the  behavior  of  a  laminated  composite,  it  is  essential  to  understand  the 
behavior  of  a  lamina.  Linear  elastic  behavior  is  typically  assumed  with  the  macro- 
mechanical  approach.  In  this  study,  the  nonlinear  material  behavior  of  Gr/PEEK  has  been 
modeled  as  increments  of  linear  behavior.  Therefore,  the  macromechanical  approach  is 
applicable. 

If  we  start  with  the  most  general  case  of  the  stress  strain  relations  for  an  anisotropic 
material,  we  can  make  use  of  the  generalized  expression  for  Hooke's  Law  (equation  1). 

<Ti  =  Cij  ej  y-l,-,6  (1) 

where: 

oj  =  stress  components 

Cij  =  stiffness  matrix 

ej  =  strain  components 


2-5 


The  convention  for  this  contracted  notation  is  the  same  as  that  used  by  [  12].  For 
clarity,  the  following  relationships  between  the  tensor  notation  and  the  contracted  notation 
are  given. 

TENSOR  NOTATION 

<>11 

°22 
O33 

*23 

% 

*12 
ell 
*22 

Y23  =  2e23 

Y3i  =  2e3i 
^12*^12 

where: 

Yy  -  engineering  shearing  strain 
£jj  =  tensor  shearing  strain 

The  linear  strain-displacement  relationships  in  contracted  notation  can  be  formally 
defined  as: 
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where:  u,v,  and  w  are  x,y,  and  z  displacements 
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The  stiffness  matrix  (Cjj)  has  36  constants.  These  36  constants  are  not  independent  and 
for  elastic  materials  the  stiffness  matrix  is  found  to  be  symmetric.  Therefore,  only  21  of 
the  constants  are  actually  independent.  The  stress- strain  relations  for  an  anisotropic  material 
(21  independent  constants)  can  therefore  be  expressed  as: 
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Equation  (2)  is  the  most  general  form  of  the  stress-strain  relations  of  an  anisotropic 
material  with  no  planes  of  material  symmetry.  Equation  (2)  can  be  simplified  for  materials 
which  have  symmetry  of  material  properties.  For  orthotropic  materials  there  are  three 
orthogonal  planes  of  material  symmetry.  Therefore,  equation  (2)  can  be  simplified  as 
shown  in  equation  (3).  For  the  orthotropic  case,  there  are  nine  independent  constants,  and 
there  is  no  interaction  between  normal  stresses  and  shearing  strains  and  vice  versa. 
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If  we  choose  to  invert  equation  (1)  and  write  an  equivalent  strain-stress  relationship,  we 
must  use  the  compliance  matrix  (Sjj)  as  follows: 

£j  -  Sjjaj  ij  =  1,  — ,6  (4) 

For  most  test  purposes,  the  compliance  matrix  components  (Sjj)  are  more  easily 
determined  than  the  components  of  the  stiffness  matrix  (Cjj).  This  is  due  to  the  fact  that 

simple  tests  are  performed  with  a  single  known  load  condition  [12].  Therefore,  it  is 
common  to  write  equation  (3)  in  terms  of  the  compliance  matrix  (Sjj)  as  shown  in 

equation  (5). 
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For  an  orthotropic  material,  the  compliance  matrix  (Sy)  can  be  written  in  terms  of 
engineering  constants.  These  constants  are  obtained  from  uniaxial  tension  and  comp¬ 
ressions  tests  as  well  as  shear  tests.  These  engineering  constants  include  Young's  moduli, 
Poisson's  ratios  and  the  shear  moduli.  Therefore,  for  an  orthotropic  material,  the 
compliance  matrix  can  be  written  as  shown  in  equation  (6). 

Since  the  specimens  dealt  with  in  this  thesis  were  essentially  thin  plates,  a  state  of  plane 
stress  was  assumed.  In  order  to  make  this  assumption,  one  assumes  that  through  the 
thickness  normal  stresses  may  be  neglected.  Without  such  an  assumption,  a  three- 
dimensional  finite  element  analysis  would  be  required.  This  study  attempted  to  accu¬ 
rately  model  failure  of  Gr/PEEK  laminates  with  a  two-dimensional  model.  It  is  important  to 
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remember  that  for  certain  ply  lay-ups,  normal  stresses  may  be  significant  and  therefore  can 
not  be  neglected.  For  a  lamina  oriented  as  shown  in  Figure  2-4,  a  state  of  plane  stress  is 
defined  by  the  following  relations. 

O3  =  0  T23  =  ^  ^31  =  0 


Figure  2-4.  Unidirectional  Lamina 
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Therefore,  equation  (5)  can  be  further  simplified  as  shown  in  equation  (7). 


'ejl  ["  Sl1  S!2  0  1  r  o1 ' 
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To  obtain  the  stress-strain  relations  from  equation  (7),  one  need  only  invert  the 
compliance  matrix  (Sy),  which  results  in  the  reduced  stiffness  matrix  (Qjj)  as  shown  in 


equation  (8). 
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The  components  of  the  reduced  stiffness  matrix  (Qjj)  may  be  defined  in  terms  of  the 
components  of  the  compliance  matrix  (Sjj)  or  in  terms  of  engineering  constants  as  [12]: 


n  S22 _ E1 

~  SnS^n  i-v12v21 


Q12-- 


v12E2  _  V21E1 


S11S22‘S  12  1-V12V2I  1-V12V21 


022  = 


SllS22-S^2  1-VI2V21 


V,-  V,,  1 
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For  an  orthotropic  lamina  subjected  to  in-plane  forces,  equations  (7)  thru  (9)  form  the 
foundation  for  stiffness  and  stress  calculations. 

Typically,  in  most  applications,  the  principal  directions  of  orthotropy  do  not  coincide 
with  the  coordinate  axes  (see  Figure  2-5).  In  these  instances,  we  must  use  a  transformation 
relationship  between  stresses  and  strains  in  the  principal  material  directions  and  those  in  the 
coordinate  directions. 


Figure  2-5.  Coordinate  &  Principal  Material  Axes 
From  simple  geometric  considerations,  a  transformation  equation  can  easily  be  written  to 
express  stresses  in  an  X-Y  coordinate  system  in  terms  of  stresses  in  the  principal  material 
directions.  This  transformation  can  be  expressed  as  follows: 


'  V 

~  m2  n2  -2nm  “ 

'V 

.  - 

n2  m2  2nm 

< 

<h  > 

_  nm  -nm  m2-  n2  _ 

'  *12 ' 

(10) 
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where: 

m  =  cos  0 
n  =  sin  0 

If  we  implement  the  transformations,  we  can  generalize  equation  (8)  to  account  for  a 
lamina  in  which  the  fibers  are  not  aligned  with  the  coordinate  axes.  The  stress-strain 
expression  can  then  be  written  as  shown  in  equation  (1 1). 


a,  r= 


> 


Oil  Ql2  0l6 

Q12  Qz2  Qas 

Ql6  Q26  0fi6  J 


(ID 


The  matrix  (Qjj)  is  referred  to  as  the  transformed  reduced  stiffness  matrix.  The 
components  of  the  transformed  reduced  stiffness  matrix  are  defined  as  follows: 

0  11  =  Quin4  +  2(Qi2  +2Q66>nW  -K&n4 

QI6*<QH' <3,2 * (Q12 - 02, *  V!" 

^  26  ”  <Q.r  0,2  •2Q«)"3"' + <Q12  •  Q,,  *  2QJ™3 

0  55  =  (Qn  +Q22 * 2Qi2  _2Q^)n^m^  +Qgg(n4  +  m4) 
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where:  m  =  cos  0 


n  =  sin  0 

In  the  general  form  of  the  stress-strain  relations,  as  shown  in  equation  (1 1),  it  is 
important  to  note  that  there  is  coupling  between  shear  strain  and  normal  stresses  and  also 
between  shear  stress  and  normal  strains.  The  mechanical  behavior  of  a  lamina  can  be 
described  by  equation  (1 1).  With  this  information,  we  can  now  proceed  to  study  the 
behavior  of  a  laminated  fiber-reinforced  composite  laminate. 

As  previously  discussed,  interlaminar  shear  stresses  develop  in  composite  laminates  due 
to  the  tendency  of  each  ply  to  deform  according  to  its  respective  material  properties.  These 
interlaminar  shear  stresses  often  result  in  the  complicated  phenomenon  of  delamination. 
Typically  interlaminar  normal  stress  (c2,  Figure  2-6)  is  the  most  predominant  stress 

affecting  delamination.  Although  the  stress  analysis  used  in  this  study  does  not  take  into 
account  the  stacking  sequence  of  a  laminate,  work  has  been  done  by  Sandhu  [19,24]  in 
formulating  "Delamination  Moment  Coefficients"  which  describe  the  tendency  of  a 
composite  to  delaminate. 

If  we  generalize  equation  (1 1)  for  each  ply  in  a  laminate,  we  can  rewrite  the  expression 
as  shown  in  equation  (13).  The  behavior  of  a  given  laminate  can  be  determined  by 
summing  the  individual  behavior  of  each  lamina.  This  approach  essentially  ignores 
delamination  effects,  and  an  equivalent  stiffness  matrix  is  generated  by  the  finite  element 
program  discussed  in  Section  II B. 


Oil  $12  0l6 

["1 

p 

►  » 

0l2  Q22  Q26 

►  (13) 

k 

-  0l6  0*  0fi6  - 

L  V  > 

r*y 

k 

where:  k  =  the  k^1  ply 
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Figure  2-6.  Interlaminar  Stresses 


B.  Linear  Finite  Element  Theory. 

The  finite  element  method  is  a  numerical  scheme  used  for  solving  problems  in  continuum 
mechanics.  As  such,  its  success  depends  totally  upon  the  reliability  of  the  computer 
program  used.  To  minimize  computer  time,  a  linear  finite  element  program  (PLSTR)  was 
used  for  several  important  aspects  of  this  thesis.  This  program  was  written  by  Dr.  R.S. 
Sandhu,  and  was  used  for  modeling,  convergence  studies,  and  contour  plotting.  The 


program  uses  a  ply-by-ply  analysis  technique  and  also  takes  into  account  the  number  and 
orientation  of  orthotropic  plies  of  a  given  laminate.  Since  PLSTR  makes  use  of  constant 
strain  triangle  elements,  the  theory  of  finite  elements  as  it  applies  to  this  element  will  be 
briefly  discussed. 

The  constant  strain  triangle  has  six  degrees  of  freedom  as  shown  in  Figure  2-7.  This 
element  is  one  of  the  earliest  finite  elements  and  is  easy  to  formulate  [6,14].  For  the 
constant  strain  triangle,  the  assumed  displacement  field  is  linear  in  the  x  and  y  coordinates. 
Using  the  convention  of  [6],  the  displacement  fields  may  be  written  as: 

u  =  [X]{at  a2  a3]  and  v  =  [X]{a4  a5  (14) 


where: 

[X]  =  [1  x  y] 


y.v 


Figure  2-7.  Constant  Strain  Triangle 
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If  we  consider  equation  (14)  at  the  three  nodes  of  Figure  2-7,  we  can  write,  for  the  u 
displacement  field,  the  following  equation. 


rur 

ran 

*  “2 

>  =  [  A]  < 

a,  l 
32  , 

.U3. 

,aj 

where: 


-  1 

X1 

yi  “ 

[A]  = 

1 

x2 

y2 

1 

x3 

y3  . 

(15) 


If  we  invert  matrix  [A]  we  obtain  the  following: 


lAI"‘=5C 


x2ft 'x3y2 
y2  -y3 
x3  -x2 


*3yt  -  *iy3 

y3  -yi 

Xj  -  x3 


xiy2  -  *2yi 

ym 

x2  -X1 


(16) 


where: 

A  =  the  area  of  the  constant  strain  triangle 
In  terms  of  nodal  degrees  of  freedom,  equation  (14)  can  be  written  as: 
u  =  [  Ni  N2  N3]  {  Ul  } 


where: 


[Nj  N2  N3]=[l  x  y  ]  [ A] 


-1 


(17) 


If  the  displacement  functions  are  differentiated,  we  obtain  the  strains  for  the  element 
(small  strain  assumption)  as  follows: 


„  du  „  dv  „  3u  3v 

^“37  Yxy*3y  +  35T 


(18) 
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The  strain-displacement  matrix  [B]  operates  on  the  displacement  matrix  {d}  to  give  us 


strains  as  follows: 


'Up 

'  ' 

<  *y 

-  =  ( B  ]  < 

U3  ► 
vi 

„  Yxy  - 

v2 

<  v3> 

where: 


t  B  ]  = 


0  1  0  0  0  0 

0  0  0  0  0  1 

L  0  0  1  0  1  0  J 


[A]'1 


(19) 


The  element  vuffness  matrix  for  an  isotropic  material  may  now  be  formulated  as  shown 
in  equation  (20). 

[k]=JA  [b]T[e]  [B]tda  =  [B]T[e]  [B]tA  (20) 


where: 


t  =  the  element  thickness 

1  v  0 


[E]« 


1-v2 


v  1  0 


(plane  stress  assumption) 


0  0 


¥  J 


Since  we  have  already  developed  an  expression  for  the  stress-strain  relation  for  an 
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orthotropic  material  (see  equation  (13)),  we  can  say  that  [Q  jj]  is  analogous  to  [E]. 
Therefore,  we  can  rewrite  equation  (20)  as  follows: 


M  =  JA  [B]T  [Q]  [B]  t  dA  =  [B]T  (Q1  [B]  t  A 


where: 

[Q]  = 


Oil  0l2  0l6 
Ql2  O22  Q26 
O16  Q*  066  - 


(21) 


The  program  PLSTR  makes  use  of  four  constant  strain  triangles  to  produce  the  quad¬ 
rilateral  elements  which  were  used  in  this  study  (see  Figure  2-8).  Since  there  are  two 
degrees  of  freedom  per  node,  the  resulting  stiffness  matrix  is  a  10  x  10.  However,  two  of 
the  rows  and  columns  of  this  matrix  are  based  on  the  fictitious  node  "M"  as  shown  in  Figure 
2-8.  Therefore,  the  stiffness  matrix  can  be  reduced  to  an  8  x  8  by  assuming  that  there  are  no 
applied  loads  at  "M".  The  [B]  matrix  may  be  reduced  from  a3xiOtoa3x8by  assuming 
that  displacements  at  node  "M”  are  given  by: 

Uw-  m+ujyK+m. 
vMa  vi+Vjy.K+vL 

In  order  to  account  for  through  the  thickness  ply  orientations,  an  equivalent  stiffness 
matrix  was  formulated.  The  program  PLSTR  makes  use  of  tite  assumption  that  dis¬ 
placements  through  the  thickness  of  the  laminate  are  constant  Therefore,  the  equivalent 
stiffness  matrix  is  formulated  by  summing  up  the  individual  stiffness  matrices  for  each  ply 
as  shown  in  equation  (22). 
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n 

[keq]  =  AX  [B]T  [QiHBlti  (22) 

i=l 

where: 

i  =  the  i^1  lamina 

n  =  the  total  number  of  plies 
tj  =  the  thickness  of  the  i*  ply 

Since  elements  of  different  plies,  through  the  thickness  of  the  laminate,  share  common 
nodes  (see  Section  HI  B  (3)),  they  will  therefore  all  deform  by  the  same  amount  for  a  given 
load.  "Whether  or  not  an  element  fails  depends  on  its  own  material  properties  a<4}  on  its 
orientation  to  the  given  load  [13]." 

C.  Nonlinear  Constitutive  Relationships. 

The  nonlinear  material  finite  element  program  used  for  this  thesis  (PLSTREN)  was 
written  by  Dr.  R.S.  Sandhu.  This  progressive-ply-failure  program  predicts  both  damage 
initiation  and  propagation  in  composite  laminates.  The  program  is  modular  by  design  and 
allows  for  the  use  of  various  finite  elements,  ply  failure  criteria,  and  post-failure  element 
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unloading  options  [21].  Although  this  program  is  well  documented  [23,25],  the  pertinent 
theory,  as  it  applies  to  this  thesis,  will  be  discussed. 

As  previously  mentioned  in  Section  I,  Gr/PEEK  exhibits  nonlinear  behavior  (also  see 
Appendix  A).  Therefore,  it  is  appropriate  to  use  an  incremental  form  of  the  constitutive  law. 
To  define  this  incremental  form  of  the  constitutive  law,  two  assumptions  are  made  [25]. 

(1)  increments  of  strain  depend  upon  the  strain  state  and  the 
increment  of  stress 

(2)  the  increment  of  strain  is  proportional  to  the  increment  of 
stress 

Therefore,  the  incremental  constitutive  law  for  an  anisotropic  material  (plane  stress 
assumption')  is  written  as: 

dei  =  Sij(ei)daj  (ij  =  1,2,6)  (23) 


Equation  (23)  can  be  written  in  matrix  form  as  follows: 


'dEi' 

"  S11  S12  S16  “ 

do^ 

<  d62 

►  s 

S21  S22  S26 

- 

do^  ► 

-  Si6  S26  Sgg  _ 

.do6. 

Equation  (24)  can  be  written  more  simply  as: 


(24) 


where: 


dOj 

*1 


d02 

dE2 

do6 

dE6 


[de]  =  [S]  [do] 


(25) 


normal  stress  increment  in  the  liber  direction 
normal  strain  increment  in  the  fiber  direction 
normal  stress  increment  in  the  transverse  direction 
normal  strain  increment  in  the  transverse  direction 
shear  stress  increment 
shear  strain  increment 

compliance  matrix  components  which  represent  the 
average  values  during  the  increment  of  stresses 


For  an  orthotropic  lamina  (see  Section  n  A),  we  can  simplify  equation  (24)  by  setting 


S16  =  S26  =  0 


(26) 


This  assumption  requires  that  the  lamina  remain  orthotropic  at  any  given  level  of  load. 
Furthermore,  we  can  define  the  remaining  components  of  the  compliance  matrix  in  terms  of 
engineering  constants  as  follows: 


S 12  -  S21  = 


.Yu.  ?i l 

En  ^22 


(27) 


If  we  combine  equations  (24),  (26),  and  (27),  we  obtain  the  incremental  stress-strain 
relations  for  a  given  lamina. 

do, 

d£i  =  p'-*'  (l-v^R)  (28) 

Ell 


where: 

R=^2.  (doi  *  0) 
dci 


Inspection  of  equations  (28)  and  (29)  reveal  that  errors  would  result  from  using  dei  or 
d£2  to  determine  the  incremental  elastic  constants  (Ei  1  or  E22)  from  stress-strain  curves 
obtained  under  simple  load  conditions.  For  the  simple  case  of  uniaxial  loading  of  an 
anisotropic  material,  the  resulting  stress  state  is  not  the  same  as  that  obtained  under  biaxial 


2-21 


loading  conditions.  This  can  be  explained  more  clearly  in  pictorial  form  where  it  is  clear  that 
simple  loading  conditions  do  not  accurately  represent  a  state  of  biaxial  stress. 


For  example,  d£2  of  equation  (29)  corresponds  to  curve  ON 
(Figure  2-9)  on  the  plane  OEHG,  while  simple  stress  strain  curve 
OM  lies  on  the  plane  OEDC  [25]. 


Since  no  experimental  data  is  available  for  the  biaxial  state  of  stress,  we  must  somehow 
take  into  account  the  effect  of  a  biaxial  stress  state.  Sandhu  accomplishes  this  by  assuming 
that  simple  equivalent  strain  increments  can  be  computed  as  follows. 


del 


6(1  1-vi2R 


(31) 


del 


(32) 


If  we  invert  equation  (23)  and  write  an  equivalent  expression  for  the  kth  ply  we  obtain 
equation  (33). 

[do]k  =  [C]k  [de]k  (33) 

where: 

[C]k  =  stiffness  matrix  of  the  ply 

[da]k  =  stress  increment  in  the  k^1  ply  relative  to 
the  material  axes  1,2 

[d£]k  =  strain  increment  in  the  k^1  ply  relative  to 
the  material  axes  1,2 


As  discussed  in  Section  n  A.  and  shown  graphically  in  Figure  2-5,  the  coordinate  axes 
(X  &  Y)  do  not  in  general  coincide  with  the  principal  material  axes.  Therefore,  the  stress 
and  strain  increments  in  the  two  coordinate  systems  must  be  related  by  an  appropriate 
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Figure  2-9.  £  2  Under  a  Biaxial  State  of  Stress  [23,25] 


transformation  matrix.  Following  the  convention  used  by  [23,25],  we  can  relate  the  stress 
and  strain  increments  in  the  two  coordinate  systems  by  the  transformation  matrix  [T]k. 


[da]k  =  [T]k  [da]k 

[d£]k  =  [T]k[de]k 


(34) 


where: 


[do]k  =  stress  increment  in  x,y  coordinate  system 
[de)k  -  strain  increment  in  x,y  coordinate  system 


Substituting  equation  (34)  into  equation  (33)  yields: 


[T]k  [da]k  =  [C]k  [de]k  (35) 

where: 

[C]k  =  [T]'1[C]k[T]k 


Since  the  assumption  is  made  that  stresses  are  uniformly  distributed  through  the 
thickness  of  each  ply,  we  can  define  the  stress  resultant  increments  [dN]  (in  the  x,y 
coordinate  system)  as  follows: 


P 


[dN]  =  X  [do]ktk 

k=l 


(36) 


where: 

tk  =  thickness  of  the  kth  ply 
P  =  total  number  of  plies  in  the  laminate 


Making  use  of  our  transformations  (equation  35),  the  stress  resultant  increments  may  be 
written  as: 


P 


*k  [  C]k  [de]k 


(37) 


Since  the  strain  increments  [d£]k  are  the  same  for  all  plies  through  the  thickness  of  the 
laminate,  we  can  write  the  following: 


[de]k  =  [de] 


Therefore  equation  (37)  becomes: 

[dN]  =  [F]  [de] 


(38) 
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where: 


P 

[F]  =  X  ^  ^ 

k=l 

If  equation  (38)  is  inverted,  [F]'1  represents  the  average  compliance  properties  of  the 
given  laminate  during  the  (n  )&  load  increment  [25]. 

[de]  =  [F]’1  [dN]  (39) 

The  problem  in  evaluating  equation  (39)  is  that  the  properties  (values)  of  [F]"1  are  not 
known  at  the  start  of  the  (n  )*h  load  increment.  Therefore,  the  properties  of  [F]*1  at  the  end 
of  the  (n-l)^1  load  increment  are  used.  This  method  is  often  referred  to  as  the  "predictor 
corrector"  method. 

[de]„=[F]-1n.,[dN]n  (40) 

These  equations  set  up  the  incremental  scheme  under  which  PLSTREN  operates.  Once 
the  strain  increments  have  been  calculated  from  equation  (40),  [dcr]*,  [de]*,  and  del^  are 

calculated.  These  stress  and  strain  increments  are  added  to  the  n1*1  load  increment  stresses 
and  strains  to  yield  the  current  stress  and  strain  states  in  all  plies.  These  "current"  stresses 
and  strains  are  used  to  calculate  the  average  elastic  properties  of  the  plies  and  a  new  [F]*1  is 
calculated.  The  incremental  scheme  is  continued  until  the  values  of  [de]n+1  and  [de]n  have 

converged  (see  equation  41). 

[de]n+i  -  [de]n 

- — - ! —  <  0.001  (41) 

[de]n 

PLSTREN  requires  that  the  mechanical  properties  of  the  laminate  be  entered  in  the  form 
of  tabular  data.  These  properties  were  determined  by  [13]  through  simple  tests  (uniaxial 
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tension  and  compression  along  and  transverse  to  the  material  axes  and  shear).  Once  this 
basic  stress-strain  data  has  been  read  by  the  program,  piecewise  cubic  spline  interpolation 
functions  [2]  are  used  to  represent  the  data.  The  advantage  of  the  cubic  spline  functions  is 
that  they  yield  smooth  stress-strain  curves  from  which  accurate  moduli  of  elasticity  over  the 
entire  range  of  the  curves  can  be  determined  [23]. 

For  clarity,  the  basic  property  tests,  resulting  curves,  and  corresponding  basic  properties 
are  shown  in  Table  2-1.  In  addition,  the  material  property  curves  are  shown  in  Appendix  A. 
As  might  be  expected,  a  laminate  can  not  continue  to  be  loaded  without  failing.  Therefore, 
the  failure  criterion  used  by  PLSTREN  is  discussed  in  Section  II D. 


Table  2-1.  Property  Tests/Resulting  Curves/Corresponding  Basic  Properties 

TEST 

CURVE 

PROPERTY  OBTAINED 

o 

Oj  vs  £j  — 

* 

0  Tension 

V 12  vs  £i  - 

^  T 

^  V12 

0 

0  Compression 

G J  VS  e  1  - 

—  E? 

v  vs  e  — 
12  1 

^  V 12 

o 

90  Tension 

T 

&2  vs  e2 

E  2 

0 

90  Compression 

a2  e2 

C 

e2 

0 

±45  Tension 

T  12  VS  ^  12 

G  12 
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D.  Strain  Energy  Failure  Theory. 

As  expected,  continued  loading  of  a  laminate  will  cause  failure.  This  failure  condition, 
under  general  states  of  stress,  is  determined  by  a  total  strain  energy  failure  criterion 
developed  by  Sandhu  [21,22,23,25].  For  nonlinear  materials,  the  failure  state  is  assumed 
to  be  a  function  of  both  stresses  and  strains.  Therefore,  a  scalar  function,  f,  is  defined 
which  accounts  for  the  failure  condition  of  a  given  material  as: 

f(o,e,K)  =  1  (42) 

where: 

cr  =  stress  state 
e  =  strain  state 
K  =  material  characteristics 

If  equation  (42)  is  written  in  an  explicit  form,  whereby  the  scalar  strain  energy  is  used  as 
a  measure  to  determine  the  effects  of  both  stress  and  strain  states  of  the  material,  equation 
(42)  can  be  expressed  (for  an  orthotropic  material)  as: 

Kij[Je.j<Tijdeij  =1  (ij  =  1,2,3)  (43) 

where: 

Ey  =  current  strain  components 

m  =  parameter  defining  the  shape  of  the  failure  surface  in 
the  strain-energy  space 

The  failure  criterion  is  based  upon  the  total  strain  energies.  For  our  special  case  of  plane 
stress,  equation  (43)  can  be  written  as: 

k[W  l  +  K2[/e2o2de2f2  +  K6[Jr%oSie6]m6  =  l  (44) 

One  will  notice  from  equation  (44)  that  the  failure  criterion  is  based  upon  the  concept  that 
strain  energies  under  longitudinal,  transverse,  and  shear  loadings  are  independent  para- 
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meters.  Therefore,  the  ratios  of  current  energy  levels  (due  to  longitudinal,  transverse,  and 
shear  loading)  to  the  maximum  energy  levels  that  are  available  for  each,  must  be  considered. 
The  shape  of  the  failure  surface  in  strain-energy  space  is  governed  by  the  shape  factor,  mj 
[20,21].  For  m1=m2=m^=2,  the  failure  surface  is  spherical,  and  for  m1=m2=m6=l,  the 
failure  surface  is  pyramidal.  This  failure  criterion  is  compared  (for  varying  mj's)  with 

several  other  theories  in  Figure  2-10  [23].  It  should  be  pointed  out  that  this  comparison  was 
confined  to  the  first  quadrant  of  stress  space,  and  used  boron-epoxy  material  systems.  For 
purposes  of  this  thesis,  and  for  lack  of  reliable  experimental  data  for  Gr/PEEK  laminas 
under  biaxial  stress  states,  mj  was  taken  to  be  unity.  Therefore,  when  the  sum  of  the  three 

ratios  of  equation  (44)  equals  unity,  the  lamina  degrades  completely.  The  next  question 
which  must  now  be  addressed  is;  does  the  fiber  fail  or  does  the  matrix  fail,  and  how  do  we 
unload  failed  elements? 

Since  equation  (44)  indicates  when  lamina  failure  occurs,  we  now  need  only  determine 
the  failure  mode.  At  present,  two  failure  modes  are  possible.  These  include  either  matrix  or 
fiber  failure.  The  matrix  failure  mode  can  occur  without  causing  fiber  failure.  When  matrix 
failure  occurs,  the  failed  lamina  continues  to  carry  loads  parallel  to  the  fibers  but  is  unloaded 
in  transverse  tension  and  shear.  If  fiber  failure  occurs,  the  lamina  is  considered  unable  to 
carry  any  load,  and  loads  are  transferred  to  the  adjacent  lamina.  To  determine  which  of  the 
two  failure  modes  are  appropriate,  Sandhu  developed  a  criterion  given  by: 

=  Check  (45) 

i=(  1.2,3) 

If  Check  2: 0. 1  then  Fiber  Failure 
If  Check  <  0. 1  then  Matrix  Failure 
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Once  ply  failure  has  occured  in  a  given  element,  the  failed  element  can  no  longer  sustain 
loads.  Therefore,  the  load  which  was  sustained  by  the  failed  element  must  be  redistributed 
to  other  adjacent  elements  of  the  laminate.  This  study  made  use  of  a  gradual  unloading 
scheme  (option  1). 

Option  1.  (Gradual  Unloading)  -  With  this  option  it  is  assumed  that  as  elements  fail, 
adjacent  elements  are  able  to  carry  the  load  previously  carried  by  the  failed  elements.  This 
unloading  scheme  is  accomplished  by  assigning  negative  values  to  the  affected  moduli  while 
continuing  the  incremental  loading  (see  Figure  2-1 1).  When  the  stresses  in  the  failed 
elements  are  reduced  to  zero,  the  moduli  are  set  to  nominally  small  values  (100  psi). 


Figure  2-11.  Gradual  Unloading  Scheme 


m.  Analysis 


As  previously  discussed  in  Section  n,  this  study  made  use  of  both  linear  and  nonlinear 
finite  element  programs.  The  primary  purpose  of  the  linear  finite  element  analysis  was  to 
develop  a  finite  element  model  for  use  with  the  nonlinear  finite  element  program.  In 
addition,  the  linear  finite  element  analysis  was  conducted  to: 

(1)  verify  the  "convergence"  (i.e.  accuracy)  of  analytical  data  obtained  by  using 
the  finite  element  model. 

(2)  study  the  effects  of  the  eccentric  circular  discontinuity  on  the  overall  stress 
state  of  the  Gr/PEEK  laminates. 

(3)  study  boundary  condition  and  gage  length  effects  on  the  overall  stress  state  of 
the  Gr/PEEK  laminates. 

The  nonlinear  finite  element  analysis  made  use  of  the  finite  element  model  developed 
from  the  linear  analysis.  The  purpose  of  the  nonlinear  finite  element  analysis  was  to 
determine  the  initiation,  progression  of  failure,  and  failure  loads  of  the  Gr/PEEK  laminates 
studied. 

As  outlined  below,  this  section  discusses  the  analytical  portion  of  this  study. 

(A)  Specimen  Geometry 

(B)  Finite  Element  Model  Development 

(Q  Convergence  Study 

(D)  Boundary  Conditions 

(E)  Eccentricity/Boundary  Condition/Gage  Length  -  -  Effects  (Linear  Analysis) 

(F)  Nonlinear  Analysis 

A.  Specimen  Geometry. 

Before  discussing  the  development  of  the  finite  element  model,  a  brief  description  of  the 
specimens  used  during  the  experimental  portion  of  this  study  is  presented.  As  previously 
mentioned  in  Section  I,  Gr/PEEK  (16  ply)  composite  laminates  were  tested.  The  four  ply 
lay-ups  studied  were  [0°16],  I90°16],  [±45°]4s,  and  [0o/±45o/90o]2s- 

Since  this  study  involved  eccentricity  effects,  several  of  the  specimens  were  not 
geometrically  symmetric  (Figure  3-1).  Therefore,  the  finite  element  models  developed 
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were  "whole  models."  This  term,  "whole  model,"  implies  that  geometric  symmetry  could 
not  be  used  to  simplify  the  model  used  in  the  analysis.  In  some  instances  (where  the 
specimens  are  geometrically  symmetric),  it  is  possible  to  represent  the  specimens  by  finite 
element  models  in  which  only  one-half  or  one-fourth  of  the  specimen  is  modeled.  This 
simplification  (not  appropriate  in  this  study)  leads  to  finite  element  models,  with  fewer 
elements,  which  require  less  computer  time  to  analyze. 


B.  Knit 


When  constructing  a  finite  element  model,  it  is  imperitive  that  the  actual  dimensions  of 
the  specimen,  applied  loading,  and  the  physical  boundary  conditions  be  modeled  as 
accurately  as  possible.  As  discussed  in  Section  IIIA,  whole  models  were  required  in  this 
study.  Therefore,  three  separate  finite  element  models  were  required  to  represent  the  three 
different  size  specimens  shown  in  Figure  3-1.  In  addition  to  this  requirement,  glass  epoxy 
tabs  as  well  as  through  the  thickness  ply  orientations  were  incorporated  into  the  finite 
element  models.  Finally,  the  finite  element  mesh  required  refinement  in  areas  where  high 
gradients  of  stress  were  expected  [6],  This  refinement  was  required  so  that  large  jumps  in 
stress  were  avoided  as  boundaries  between  finite  elements  were  crossed.  This  was 
particularly  important  in  this  study,  since  the  quadrilateral  elements  used  for  modeling  were 
constructed  from  constant  strain  triangle  elements  (see  Section  II B). 

The  last  complication  to  the  development  of  the  finite  element  model  was  the  use  of  a 
mounting  fixture.  The  mounting  fixture  (see  Section  IV  D)  allowed  in-plane  rotation  of  the 
specimen. 

(1)  Modeling  of  Specimen  Widths. 

After  a  considerable  amount  of  work,  a  method  was  proposed  by  Dr.  Sandhu  to 
develop  one  model  which  could  be  used  to  represent  the  three  different  width  specimens 
(Figure  3-1).  This  method  involved  the  development  of  a  finite  element  model  (Figure  3-2) 


3-2 


Figure  3-1.  Specimen  Geometry 


which  represented  the  widest  specimen  shown  in  Figure  3-1. 

To  represent  the  other  two  (more  narrow)  specimens  shown  in  Figure  3-1,  the  finite 
element  model  (Figure  3-2)  was  modified.  This  modification  involved  the  mathematical 
elimination  of  the  rows  of  finite  elements  at  the  bottom  of  the  model  (Figure  3-2).  This 
was  accomplished  by  "zeroing  out"  the  elements  below  rows  AA  or  BB  (Figure  3-2).  The 
term  "zeroing  out"  an  element  implies  that  a  nominally  low  modulus  of  1000  psi  (compared 
to  the  actual  modulus  of  19,200,000  psi)  was  assigned  to  that  particular  element.  There¬ 
fore,  the  zeroed  out  elements  contributed  an  insignificant  amount  to  the  overall  stiffness  of 
the  model.  By  reducing  the  model  shown  in  Figure  3-2,  it  was  possible  (by  zeroing  out 
rows  of  elements)  to  develop  the  finite  element  models  shown  in  Figure  3-3.  The 
significance  of  this  reduction  was  that  one  finite  element  model  was  used  to  represent  three 
different  width  specimens. 

(2)  Modeling  of  Glass  Edoxv  Tabs. 

As  previously  mentioned,  it  was  necessary  to  incorporate  glass  epoxy  tabs  (Figure  3-1) 
in  the  finite  element  model.  This  was  accomplished  by  further  modifications  to  the  finite 
element  models  shown  in  Figure  3-3.  The  addition  of  tabs  to  the  model  was  accomplished 
by  adding  an  additional  finite  element  (through  the  thickness)  to  each  element  contained 
within  the  tab  area  (see  Figures  3-4  B  and  3-4  C).  Since  out-of-plane  bending  was  not 
allowed  (analytically),  the  total  thickness  of  the  upper  and  lower  tabs  (Figure  3-4  A)  were 
modeled  as  one  element  with  a  thickness  equal  to  the  sum  of  the  upper  and  lower  tab 
thicknesses  (Figure  3-4  C).  Also,  note  that  the  finite  element  mesh  is  refined  as  you  move 
into  the  tab  area  (A  to  B,  Figure  3-4  B).  This  refinement  was  incorporated  due  to  an 
expected  stress  gradient  arising  from  the  geometric  discontinuity  caused  by  the  tab  material. 
The  tabs  were  tapered  to  minimize  this  geometric  discontinuity  (Figure  3-4  A). 
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Figure  3-3.  Finite  Element  Models  (no  tabs) 
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Figure  3-4  A.  Cross-Sectional  View  of  Specimen  ( left  end) 
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Figure  3-4  B.  Finite  Element  Model 
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Figure  3-4  C.  Section  CC 

•  indicates  elements  added  to  finite  element  model  to 
account  for  the  glass  epoxy  tabs;  the  thickness  of 
these  elements  is  the  sum  of  the  upper  and  lower 
tab  thicknesses. 


Figure  3-4.  Tab  Area  Modeling 


(3)  Modeling  of  Through  the  Thickness  Plv  Orientations. 

Since  four  different  ply  lay-ups  (stacking  sequences)  were  considered  in  this  study,  it 
was  necessary  to  model  the  various  ply  orientations  through  the  thickness  of  the  laminates. 
This  was  a  simple  task  for  both  the  [O0^]  and  [90oig]  laminates.  As  shown  in  Figure  3-5, 
both  of  these  laminates  were  modeled  as  one  through  the  thickness  element.  For  the 
[±45°]4S  laminates,  the  laminate  was  modeled  as  two  through  the  thickness  elements  as 
shown  in  Figure  3-5.  Finally,  the  (0°/±45°/90°]2S  laminates  (quasi  isotropic)  were  modeled 
as  four  through  the  thickness  elements  as  shown  in  Figure  3-5. 
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Figure  3-5.  Through  the  Thickness  Finite  Element  Modeling 


Since  through  the  thickness  finite  elements  share  the  same  four  nodes  (see  Section  Q 
B,  also  Figure  3-5),  the  finite  element  models  for  the  four  stacking  sequences  each 
contained  682  nodes.  However,  since  the  four  stacking  sequences  required  a  different 
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number  of  elements  through  the  thickness  of  the  laminate  (Figure  3-5),  the  total  number  of 
elements  for  the  four  stacking  sequences  were  as  follows: 

[0oig]  831  elements 

[90°i6]  831  elements 

[±45°]4s  -  1 142  elements 

[0°/±45o/90°]2s  2664  elements 

In  order  to  facilitate  the  entry  of  required  input  data,  a  FORTRAN  program  was 
written  to  take  the  finite  element  model  developed  for  the  [0°ig]  laminates  and  make  the 
modifications  necessary  to  develop  the  models  for  the  [±45°]4s  and  [0o/±45o/90°]2s 
laminates. 

(4)  Modeling  of  the  Mounting  Fixture. 

The  last  complication  to  the  finite  element  model  was  the  use  of  a  mounting  fixture. 
Photographs  of  this  fixture  as  well  as  a  detailed  explanation  of  the  fixture  are  provided  in 
Section  IV  D.  As  shown  in  Figure  3-6  (also  in  more  detail  in  Figure  4-8),  the  fixture 
utilized  steel  blocks  to  rigidly  secure  the  specimens  during  testing.  These  steel  blocks 
were  allowed  to  rotate  (in-plane)  through  the  use  of  bearings  (see  Figures  3-6  and  4-8). 

Since  the  steel  blocks  shown  in  Figure  3-6  were  allowed  to  rotate  (in-plane)  with  the 
specimen,  the  tab  area  of  the  finite  element  model  was  modified  to  include  the  steel  blocks. 
This  was  accomplished  by  using  the  equivalent  material  properties  for  a  combination  of  the 
glass  epoxy  tabs  and  the  steel  blocks  (see  Figure  3-7).  These  equivalent  material  properties 
were  obtained  by  using  the  Program  SQ5  [18]. 

Therefore  as  shown  in  Figure  3-7,  the  steel  blocks  of  the  mounting  fixture  and  the 
glass  epoxy  tabs  (non*tapered  portion)  were  combined  to  form  through  the  thickness 
finite  element  with  the  equivalent  material  properties  of  steel  and  glass  epoxy. 
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Figure  3-6.  Specimen  Mounting  (Fixture) 
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The  next  question  which  must  be  addressed  is;  does  the  model  accurately  represent  the 
specimens  used  during  experimentation? 

C.  Convergence  Study. 

In  order  to  verify  the  suitability  of  the  finite  element  model,  the  zero  eccentricity  model 
(Figure  3-3)  was  used  in  a  convergence  study.  The  zero  eccentricity  model  was  used 
because  a  detailed  convergence  study  was  completed  by  [10]  on  a  very  similar  model,  and 
in  addition,  abundant  stress  intensity  factor  information  is  available  for  symmetric  isotropic 
plates  containing  circular  discontinuities  [15].  Therefore,  the  zero  eccentricity  model  was 
given  isotropic  properties  as  follows: 

Ei  i  =  E22  =  19,200,000  psi 


v  =  .3 
Gi2  =  ; 


ill 


■=  7,384,615  psi 


2(l+v) 

The  isotropic  model  was  then  run  using  the  linear  finite  element  program  (see  Section 
II B),  and  the  stress  intensity  factor  (SIF)  at  the  edge  of  the  hole  (point  A,  Figure  3-8  A) 
was  found  to  be  3.29.  Peterson  [15]  reports  the  SIF  as  3.45,  and  Fisher  [10]  found  the 
SIF  to  be  3.28.  The  goal  of  this  convergence  study  was  to  insure  that  the  finite  element 
model  would  give  a  SIF  within  5%  of  published  values  for  the  theoretical  value  of  the  SIF. 
A  comparison  of  the  SIF  obtained  from  this  model  and  the  SIF  reported  by  Peterson  [15], 
showed  a  percent  error  of  4.6%.  From  this  comparison  of  SIFs,  and  from  the  mesh 
refinement  studies  completed  by  [10],  the  mesh  refinement  of  the  finite  element  model  was 
found  to  be  satisfactory. 


D.  Boundary  Conditions- 

Once  the  finite  element  models  had  been  developed  (Section  HI  B),  the  appropriate 
boundary  conditions  were  required  to  complete  the  modeling  process.  Since  a  mounting 
fixture  (Section  m  B(4))  was  used  in  this  study,  two  sets  of  boundary  conditions  were 
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required.  First,  boundary  conditions  were  developed  for  models  which  did  not  allow  in¬ 
plane  rotation.  And  secondly,  a  set  of  boundary  conditions  were  developed  for  the  models 
which  allowed  in-plane  rotation. 

(1)  Boundary  Conditions  (No  Mounting  Fixture!. 

As  shown  in  Figure  3-8  A,  boundary  conditions  were  developed  for  the  models  which 
did  not  allow  in-plane  rotation.  These  boundary  conditions  were  as  follows: 

a) .  Nodes  located  along  the  line  x=3.8"  were  given  prescribed 

displacements  of  zero  in  both  the  x  and  y-directions. 

b) .  Nodes  located  along  the  line  x=10.2"  were  given  a  prescribed  displace¬ 

ment  in  the  x-direction,  and  a  zero  displacement  in  the  y-direction. 

c) .  All  other  nodes  in  the  model  were  free  to  displace  in  both  the  x  and  y- 

directions. 

These  boundary  conditions  were  developed  to  simulate  the  actual  boundary  conditions 
present  during  experimentation  without  the  mounting  fixture.  These  boundary  conditions 
resulted  from  specimen  mounting  in  the  grips  of  the  Instron  Test  Machine  (Section  IV  C, 
Figure  4-6).  It  is  important  to  no.e  that  prescribed  displacements  were  not  given  to  nodes 
which  surrounded  "zeroed  out"  elements  (see  Section  in  B(l)). 

(2)  Boundary  Conditions  (Mounting  Fixture). 

As  shown  in  Figure  3-8  B,  boundary  conditions  were  developed  for  models  which 
allowed  in-plane  rotation.  These  boundary  conditions  were  as  follows: 

a) .  The  node  (node  I)  located  at  the  geometric  centerline  of  the  model,  and  at 

x=3.8",  was  given  prescribed  displacements  of  zero  in  both  the  x 
and  y-directions. 

b) .  The  node  (node  J)  located  at  the  geometric  centerline  of  the  model,  and  at 

x*  10.2",  was  given  a  prescribed  displacement  in  the  x-direction,  and 
a  zero  displacement  in  the  y-direction. 

c) .  All  other  nodes  in  the  model  were  free  to  displace  in  both  the  x  and  y- 

directions. 

These  boundary  conditions  were  developed  to  simulate  the  actual  boundary  conditions 
present  during  experimentation  with  the  mounting  fixture.  Since  in-plane  rotation  of  the 
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Figure  3-8  B.  Boundary  Conditions  (mounting  fixture) 


Point  "o"  is  the  origin  of  the  coordinate  system. 


Figure  3-8.  Boundary  Conditions  (Modeling) 
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specimen  was  allowed,  only  two  nodes  (nodes  I  and  J,  Figure  3-8  B)  were  given 
prescribed  displacements.  Therefore,  in-plane  rotation  was  allowed  about  axes  (normal  to 
the  plane  of  Figure  3-8  B)  through  nodes !  and  J.  For  a  detailed  explanation  of  the 
mounting  fixture,  as  well  as  photographs,  refer  to  Section  IV  D.  It  is  important  to  note  that 
as  the  specimen  widths  changed  (Figure  3-1),  the  y-location  of  the  geometric  centerline  also 
changed. 

E.  Eccentricitv/Boundarv  Condition/Gage  Length  -  -  EfffiflS  (Linear  Analysis)- 

As  previously  discussed,  the  linear  analysis  involved  more  than  just  the  development  of 
a  suitable  finite  element  model.  An  extensive  linear  study  was  completed  before  starting  the 
nonlinear  finite  element  analysis  (Section  HI  F).  This  study  was  undertaken  to  gain  insight 
into  eccentricity,  boundary  condition,  and  gage  length  effects.  As  outlined  below,  this 

section  discusses  the  results  from  the  linear  analysis. 

(1)  Contour  Plotting 

(2)  Graphs  (Stress  States) 

a)  Near  Held/Far  Field  Stress  States 

b)  Uniformity  "  T'"r  Field  Axial  Stress  (gage  length  effects) 

(1)  CpntQUr  PiQl.ti.ng. 

To  study  the  effects  of  both  the  eccentric  circular  discontinuity  and  the  mounting  fixture 
on  the  overall  stress  states  of  the  Gr/PEEK  laminates,  contours  were  plotted.  The  output 
strain  data  obtained  from  the  linear  finite  element  program  was  used  to  plot  contours  for 
the  specimens  tested  experimentally  (see  Tables  4-1  through  4-4).  These  contour  plots 
(Appendix  B)  are  grouped  according  to  the  stacking  sequence  of  the  laminate,  and  also  by 
the  boundary  conditions  used  with  the  finite  element  models  (see  Section  in  D).  In  order 
to  discuss  the  results  from  this  study,  the  contour  plots  for  the  [90°i5]  laminates  have  been 
taken  from  Appendix  B  and  are  shown  in  Figures  3-9  and  3-10.  The  terms  "hinge"  and 
"no  hinge"  in  these  figures  indicates  that  in-plane  rotation  (analytically)  was  or  was  not 
allowed,  respectively  (see  Section  HI  D). 
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Rgore3-9.  Contour  Plots  (90  degree  lay-up  /  no  mounting  fixture) 
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Figure  3-10.  Contour  Ptott  (90  degree  Uy-up  /  mounting  fixture) 
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Close  inspection  of  Figures  3-9  and  3-10  provide  important  information.  First,  note 
that  the  contours  are  not  symmetric  for  the  nonsymmetric  models  (eccentricity  *  0.0").  The 
contours  show  higher  levels  of  strain  above  the  discontinuity  compared  with  the  levels  of 
strain  (at  an  equal  distance)  below  the  discontinuity.  This  verified  an  early  assumption  that 
the  eccentric  circular  discontinuity  would  create  a  centroid  (area)  shift,  and  therefore  the 
tensile  loading  (from  the  Instron)  would  produce  a  bending  moment  (see  Figure  3-11). 

This  bending  moment  (Figure  3-11),  produced  a  positive  (tensile)  stress  above  the 
discontinuity  and  a  negative  (compressive)  stress  below  the  discontinuity.  The  stresses 
from  the  bending  moment,  combined  with  the  tensile  loading  from  the  Instron,  produced 
higher  strains  (for  the  nonsymmetric  models)  above  the  discontinuity.  Secondly,  no 
apparent  differences  were  noted  between  Figures  3-9  and  3-10.  Even  if  Figures  3-9  and  3- 
10  were  placed  on  top  of  each  other  (overlaid),  only  very  slight  differences  between  the 
contours  were  noted.  This  is  significant  because  Figure  3-9  shows  contours  for  finite 
element  models  which  do  not  allow  in-plane  rotation;  while  Figure  3-10  shows  contours 
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(centroid  for  nonsymmetric  model)  As  =  Centroid  (area)  Shift 


The  eccentric  circular  discontinuity  causes  a  centroid  (area)  shift  Therefore, 
the  applied  load  (from  the  Instron)  creates  a  bending  moment  (m  )• 


Mf  =  A  S  X  Load 


Figure  3-11.  Centroid  Shift  and  Resulting  Bending  Moment 
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for  finite  element  models  in  which  in-plane  rotation  is  allowed  (see  Section  IEt  D). 

Therefore,  the  use  of  the  mounting  fixture  (see  Section  IV  D),  which  allowed  in-plane 
rotation,  did  not  appreciably  affect  the  stress  state  near  the  discontinuity. 

Finally,  contours  were  also  plotted  to  verify  the  validity  of  the  output  data  obtained  from 
the  linear  finite  element  program.  These  contours  are  compared  to  experimentally  obtained 
photographs  of  the  photoelastic  stress  states  of  the  Gr/PEEK  laminates  tested  (see  Section 
V  A).  Since  experimentation  was  conducted  for  this  comparison,  a  discussion  of  the 
results  of  this  study  are  discussed  in  Section  V  A. 

(2)  Graphs  (Stress  States’), 
a)  Near  Field/Far  Field  Stress  States. 

In  addition  to  contour  plots,  graphs  of  both  near  field  (Figure  3-12)  and  far  field  (Figure 
3-13)  stress  states  were  constructed  from  the  output  strain  data.  These  graphs  were 
constructed  to  study  boundary  condition  and  eccentricity  effects  on  both  the  near  field  and 
far  field  stress  states.  Strains  were  used  for  these  graphs  instead  of  stresses  because  strains 
were  continuous  through  the  thickness  of  the  laminates  (see  Section  II B).  These  graphs 
(Appendix  C)  are  grouped  according  to  the  sucking  sequence  of  the  laminate,  and  also  by 
the  boundary  conditions  used  with  the  finite  element  models  (see  Section  in  D).  In  order 
to  discuss  the  results  from  this  study,  the  graphs  for  the  [90oig]  laminates  have  been  taken 
from  Appendix  C,  and  are  shown  in  Figures  3-14  and  3-15. 

Inspection  of  Figures  3-14  and  3-15  provide  important  information.  First,  Figure  3-14 
shows  the  near  field  axial  strains  for  the  [90°i5]  laminates.  Note  that  the  strains  (for  the 
nonsymmetric  models)  are  not  symmetric,  with  the  strain  above  the  discontinuity  being 
higher  than  the  strain  (at  an  equal  distance)  below  the  discontinuity.  This  can  be  attributed 
to  the  bending  moment  (Figure  3-1 1)  discussed  in  Section  III  E  (1).  Also,  note  that  there  is 
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Figure3-12.  Near  Field  Location  for  Figure  3-14 
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Figure  3-13.  Far  Field  Location  for  Figure  15 
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no  apparent  difference  between  the  graphs  (Figure  3-14)  showing  data  with  and  without  the 
mounting  fixture.  This  again  suggests  that  the  use  of  the  mounting  fixture  (Section  IV  D) 
did  not  affect  the  stress  state  near  the  discontinuity.  Finally,  Figure  3- IS  shows  the  far 
field  axial  strains  for  the  [90°16]  laminates.  These  graphs  (Figure  3-15)  show  that  the 
variations  in  far  field  axial  strains  were  reduced  by  the  use  of  mounting  fixture.  However, 
as  already  noted  in  Figure  3-14,  this  reduction  in  the  variations  of  the  far  field  axial  strains 
did  not  affect  the  state  of  stress  near  the  discontinuity.  Therefore,  although  the  use  of  the 
mounting  fixture  caused  a  more  uniform  state  of  far  field  stress,  the  linear  analysis 
predicted  that  the  near  field  stress  state  would  not  be  affected  by  the  use  of  the  mounting 
fixture. 

b)  Uniformity  of  Far  Field  Axial  Stress  (gage  length  effects!. 

References  [10,13]  had  assumed,  in  their  studies  of  Gr/PEEK,  that  a  uniform  state  of 
stress  existed  far  field  (see  Figure  3-13).  Therefore,  a  study  was  conducted  to  determine  if 
in  fact  a  uniform  state  of  stress  existed  far  field.  Since  the  effects  of  the  discontinuity 
would  be  reduced  at  the  far  field  location  as  the  gage  length  (length  between  the  tabs)  of  the 
specimens  was  increased,  a  model  with  a  12"  gage  length  was  developed  (Figure  3-16). 
With  this  model,  the  variation  in  axial  stress  (for  a  constant  value  of  x)  was  studied  (see 
Figure  3-16). 

Figure  3-17  shows  the  results  from  this  study.  Note  that  the  far  field  axial  stress  is 
uniform  at  x=0"  (see  Figure  3-17).  This  corresponds  to  a  gage  length  of  12"  as  shown  in 
Figure  3-16.  As  the  discontinuity  was  approached  (x  increased.  Figure  3-16)  the  variations 
of  axial  stress  became  larger  (Figure  3-17).  Therefore,  a  12"  gage  length  would  be 
required  to  insure  that  a  uniform  state  of  far  field  stress  exists.  This  was  a  problem 
because  the  Gr/PEEK  panels  available  for  this  thesis  were  only  large  enough  to  fabricate 
specimens  with  6"  gage  lengths  (same  size  as  used  by  [10,13] ).  Therefore,  since  failure 


Note.  This  model  was  used  to  generate  the  data  shown  in  Figure 
3-17.  The  variations  in  axial  stress  along  x=constant  lines  (such 
as  A,  B,  and  C)  are  plotted  in  Figure  3-17. 


Figure  3-16.  12"  Gage  Length  Model 
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See  Figure  3-16  for  the  meaning  of  "x  distance" 


Figure  3-17.  Variation  of  Axial  Stress  versus  Distance 


would  initiate  at  the  discontinuity,  a  study  was  conducted  to  see  if  gage  length 
effects  would  affect  the  stress  intensity  factor  (SIF)  at  the  discontinuity. 

To  see  if  the  SIF  at  the  discontinuity  would  be  affected  by  the  gage  length,  three 
different  models  were  constructed  (see  Figure  3-18).  These  models  had  gage  lengths  of 
12",  10",  and  6".  Linear  computer  runs  were  made  for  these  models,  and  the  SIFs  at  the 
discontinuity  (point  A,  Figure  3-18)  were  calculated.  For  all  three  models,  the  SIFs  (at 
point  A)  were  found  to  be  equal  to  three  significant  digits  (SIF=3.29).  Therefore,  the  gage 
lengths  did  not  affect  the  state  of  stress  near  the  discontinuity.  As  a  result,  a  6"  gage 
length  was  also  used  for  this  thesis. 


Figure  3-18.  Three  Different  Gage  Length  Models 


F.  Nonlinear  Analysis. 

As  previously  mentioned,  Gr/PEEK  exhibits  nonlinear  stress-strain  behavior  (see 
Appendix  A).  Therefore,  to  accurately  study  Gr/PEEK  laminates,  a  nonlinear  analysis 
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was  required.  The  nonlinear  material  finite  element  program  used  for  this  thesis  (discussed 
in  Section  II C)  requires  the  use  of  efficient  computer  resources. 

To  make  a  nonlinear  run,  a  batch  file  was  used  to  route  the  program,  and  the  required 
input  files,  from  a  CYBER  mainframe  computer  to  a  CRAY  mainframe  computer.  The 
program  was  then  executed  on  the  CRAY.  Past  experience  with  this  nonlinear  program 
[10,13]  has  shown  that  it  is  actually  less  expensive  to  execute  the  program  on  the  CRAY 
than  it  is  on  the  CYBER.  Although  the  CRAY  performed  the  executions  at  a  higher 
expense  per  second  than  the  CYBER,  the  execution  time  was  significantly  faster,  thus 
saving  on  the  overall  expense  of  the  analysis  [13]. 

Once  a  nonlinear  run  had  been  completed  by  the  CRAY,  the  batch  file  then  routed  the 
output  files  back  to  the  CYBER.  These  output  files  are  discussed  in  Appendix  D. 

In  order  to  verify  the  nonlinear  finite  element  analysis,  experimentation  was  conducted. 
The  following  Section  (Section  IV)  outlines  the  details  of  the  experimental  portion  of  this 
study. 


iv.  Experimentation 


The  purpose  of  the  experimental  portion  of  this  study  was  to  determine  the  ultimate 
strength,  failure  modes,  and  stress-strain  response  of  Gr/PEEK  laminates  containing  an 
eccentrically  located,  0.4  in  diameter  circular  discontinuity.  Two  sets  of  experiments  were 
conducted.  In  one  set,  the  purpose  was  to  verify  the  predictions  of  the  nonlinear  finite 
element  analysis.  The  other  set  of  experiments  was  conducted  to  verify  the  validity  of  the 
analytic  stress  state  predictions  of  the  linear  finite  element  program.  The  laminates  tested 
experimentally  consisted  of  four  different  ply  lay-ups  as  follows:  [0°16],  [90°  16],  [±45°]4S, 
and  [0°/t45o/90°]2s. 

Experimentation  was  conducted  at  the  Structures  Division  of  the  Flight  Dynamics 
Laboratory  (Wright-Patterson  AFB,  Ohio).  A  total  of  fifty-two  specimens  were  fabricated. 
Forty-eight  of  these  specimens  were  instrumented  with  strain  gages  while  the  remaining 
four  were  used  for  photoelastic  testing.  All  testing  was  accomplished  at  room  temperature 
using  a  20-kip  Instron  Universal  Test  Machine  (referred  to  as  an  "Instron"). 

The  specific  objectives,  of  the  experimental  portion  of  this  study,  were  as  follows: 

(1)  Determine  the  ultimate  strength,  the  "gross"  failure  modes,  stress-strain 
response,  and  the  boundary  condition  effects  of  Gr/PEEK  laminates  containing  a  0.4  in 
diameter  (eccentrically  located)  circular  discontinuity. 

(2)  Determine  the  overall  stress  state  of  Gr/PEEK  laminates,  containing  a  0.4 
in  diameter  circular  discontinuity,  through  photoelastic  techniques. 

A.  Specimen  Description. 

The  specimens  used  in  this  study  were  fabricated  from  panels  of  Gr/PEEK  supplied  by 
the  Fiberite  Corporation  of  Great  Britain.  These  panels  consisted  of  four  different  ply  lay¬ 
ups  (previously  mentioned)  of  APC-2  with  Hercules  AS4  graphite  fibers  contained  in  a 
Victrex®  PEEK  matrix.  Fiber  content  was  guaranteed  by  the  manufacturer  to  be  61%  by 
volume  and  68%  by  weight 

All  panels  were  subjected  to  a  C-scan  evaluation  by  [10]  to  insure  that  they  were  free  of 
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significant  defects.  The  C-scan  testing  was  conducted  by  the  Non-Destructive  Branch  of 
the  Materials  Laboratory  (Wright-Patterson  APB,  Ohio). 

From  the  panels  of  Gr/PEEK,  fifty-two  specimens  were  fabricated.  Tabs  were  bonded 
to  these  specimens  as  shown  in  Figure  4-1.  The  tabs  were  constructed  from  G-10  glass 
epoxy  [0°/90°],  and  allowed  for  a  uniform  transfer  of  load  from  the  Instron  to  the  Gr/PEEK 
specimens.  In  addition,  the  tabs  provided  a  contact  area,  for  the  knurled  mounting  grips  of 
the  Instron,  to  secure  the  specimens;  this  contact  area  prevented  damage  to  the  specimens. 


The  values  of  eccentricity  used  for  this  study  (see  Figure  4-2),  were  chosen  for  two 
different  reasons.  First,  it  was  imperitive  that  the  ultimate  strength  of  the  specimens  did  not 
exceed  the  load  capacity  of  the  Instron.  And  secondly,  the  grips  utilized  by  the  Instron 
required  specimens  with  widths  less  than  two  inches.  From  Figure  4-2,  it  is  obvious  that 


Gr/PEEK 


0.397"  hole 


B.  Specimen  Instrumentation. 

As  previously  mentioned,  forty-eight  specimens  were  instrumented  with  strain  gages. 
These  gages  were  bonded  to  the  specimens  as  shown  in  Figure  4-3. 


The  two  gages  located  adjacent  to  the  discontinuity  were  referred  to  as  "near  field" 
gages,  while  the  two  gages  located  near  the  tab  area  were  referred  to  as  "far  field"  gages. 
The  location  of  the  four  strain  gages  with  respect  to  the  finite  element  model  are  shown  in 
Figure  4-4.  As  shown  in  Figure  4-4,  the  strain  gages  were  positioned  with  the  center  of 
each  gage  located  on  the  centroid  of  a  single  finite  element.  This  placement  was  used 
because  both  the  linear  and  nonlinear  finite  element  programs  give  stress  and  strain  out¬ 
put  at  the  centroid  of  the  finite  elements.  Therefore,  no  interpolation  was  required  to 
compare  experimental  and  analytical  results.  The  far  field  gages  were  used  to  determine 
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•  indicates  strain  gage  locations 

NOTE:  The  gage  near  the  tab  area  is  mirrored  by  another  strain  gage  on  the 
reverse  side  of  the  specimen  (see  Figure  4-3). 


Figure  4-4.  Strain  Gage  Locations  (Finite  Element  Model) 


the  failure  strain  for  each  laminate.  These  far  field  gages  were  placed  where  the  stresses  in 
the  laminate  were  considered  uniform.  The  analytical  study  (Section  in  E  (2)b)  showed 
that  a  far  field  "uniform  stress"  state  was  never  obtained  with  the  six  inch  gage  length 
specimens  used  during  testing.  However,  the  variation  of  far  field  axial  stresses  was 
found  to  be  small  (see  Figure  3-17).  The  near  field  gages  were  placed  as  close  to  the 
discontinuity  as  physically  possible,  while  still  insuring  that  the  center  of  the  gages  were 
located  on  the  centroid  of  a  single  finite  element 

All  four  of  the  strain  gages  were  strain  gage  rosettes  (see  Figure  4-5)  which  have  the 
capability  to  measure  longitudinal,  shear,  and  transverse  strains. 
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The  far  field  gages  were  placed  back-to-back  (i.e.  the  back  gage  mirrored  the  front 
gage).  This  configuration  insured  that  if  one  of  the  gages  failed,  one  gage  would  still 
gather  test  data  at  the  far  field  location.  In  addition,  if  both  gages  functioned  properly,  the 
far  field  strain  could  be  taken  as  the  average  values  obtained  from  both  of  the  far  field 
gages.  This  allowed  for  the  detection  of  a  misplaced  gage  or  specimen  misalignment  in  the 
Insuon.  The  near  field  gages  were  not  located  back-to-back  due  to  nonavailability  of  the 
strain  gages  during  testing. 

The  strain  gage  rosettes  were  located  so  that  the  A-legs  (see  Figure  4-5)  were  parallel  to 
the  loading  axis  of  the  specimens.  Therefore,  the  A-legs  measured  longitudinal  strain,  the 
B-legs  (in  conjunction  with  data  obtained  from  the  A  and  C  legs)  provided  shear  strain 
data,  and  the  C-legs  measured  transverse  strains.  Each  strain  gage  was  numbered  as 
shown  in  Figure  4-5.  This  numbering  scheme  was  used  for  all  of  the  specimens  so  that  no 
confusion  would  result  when  connecting  the  gages  to  the  data  acquisition  system. 

C.  Instrumented  Specimen  Testing- 

All  forty-eight  instrumented  specimens  were  tested  (in  tension)  to  their  ultimate  strength. 
The  Instron  (see  Figure  4-6)  was  used  to  apply  this  tensile  loading  under  room  temperature 
conditions.  The  loading  rate  for  all  testing  was  a  constant  crosshead  (see  Figure  4-6) 
velocity  of  0.05  inches  per  minute.  The  crosshead  is  the  component  of  the  Instron  which 
provides  the  displacement  loading  [13].  The  loading  rate  of  0.05  inches  per  minute  was 
chosen  primarily  for  the  reason  that  previous  testing  of  Gr/PEEK  indicated  that  this  rate 
was  a  "median  speed  between  a  faster  impact  type  loading  and  a  slower  creep  type  loading 
[13]." 

During  testing,  the  load  as  a  function  of  time  was  obtained  by  two  methods  as  follows: 

(1)  The  Instron  provided  strip  chart  information  (see  Figure  4-6)  of  load 
versus  time  through  the  use  of  permanently  mounted  strain  gages  which 
measured  the  crosshead  displacement 
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(2)  A  voltmeter  (see  Figure  4-6)  was  used  to  measure  resistance  across  the 
load  cell  of  the  Instron.  This  resistance  was  a  function  of  both  the  load 
setting  of  the  Instron  and  the  applied  load.  Therefore,  the  load  at  any 
time  could  be  determined  based  upon  the  reading  from  the  voltmeter. 

Since  all  testing  was  recorded  on  vMeo  tape  (see  Figure  4-6),  the  voltmeter  gave  a 
visual  indication  of  the  applied  load  at  any  time.  As  shown  in  Figure  4-6,  the  voltmeter 
was  located  adjacent  to  the  specimen.  This  setup  was  ideal  for  post-failure  analysis  since 
specimen  failure  could  be  easily  correlated  with  the  failure  load. 


A  data  acquisition  system  (see  Figure  4-7)  was  required  to  record  test  data.  This  test 
data  included  both  the  applied  load  and  the  strain  from  each  leg  of  the  four  strain  gages. 
The  load  was  determined  by  measuring  the  resistance  across  the  load  cell  during  testing. 
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For  each  leg  of  the  strain  gages,  one  pair  of  leads  (wires)  was  connected  to  a  terminal  strip 
(see  Figure  4-6).  From  the  terminal  strip,  signals  were  sent  to  the  data  acquisition  system. 
The  data  acquisition  system  consisted  of  a  strain  gage  conditioner,  amplifier,  multiplexer, 
and  VAX  computer.  The  data  acquisition  system  was  controlled  by  a  Zenith  personal 
computer. 


Figure  4-7.  Data  Acquisition  System 


Test  data  could  be  recorded  up  to  three  times  per  second.  This  rate  was  referred  to  as 
the  "sampling  rate."  Based  upon  past  studies  of  Gr/PEEK  [  13],  a  sampling  rate  of  2 
samples  per  second  was  used  for  [0°i<s],  [±45°]4S,  and  [0°/i45o/90o]2S  laminates.  For 
[90°  id  laminates,  a  sampling  rate  of  3  samples  per  second  was  used  since  these  laminates 
failed  at  relatively  low  loads.  By  increasing  the  sampling  rate  for  the  [90°i$]  laminates, 
more  data  was  collected,  and  the  strain  at  failure  was  more  accurate. 
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D.  Mounting  Fixture. 

As  previously  stated,  one  objective  of  the  experimental  portion  of  this  study  was  to 
determine  the  boundary  condition  effects  on  the  failure  modes  and  ultimate  strength  of  the 
Gr/PEEK  specimens.  Since  the  specimens  tested  in  this  study  were  not  all  geometrically 
symmetric  (see  Figure  4-2),  in-plane  bending  was  anticipated  due  to  a  centroid  (area)  shift 
created  by  the  eccentrically  located  circular  discontinuity  (see  Figure  3-11).  To  study  this 
in-plane  bending,  a  mounting  fixture  (see  Figure  4-8)  was  fabricated.  This  fixture  was 
fabricated  from  4340  steel  and  was  heat  treated  to  160  ksi  after  fabrication.  The  fixture 
design  was  based  on  a  previous  design  by  Sandhu  which  was  used  in  1985  by  Cron  [7]. 


Figure  4-8.  Mounting  Fixture 
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^  The  mounting  fixture  allowed  the  specimens  to  rotate  (in-plane).  This  rotation  was 

|  possible  because  the  specimen  was  mounted  in  steel  blocks  which  contained  bearings 

(Figure  4-9).  The  bearings  were  "seated"  on  steel  pins  which  extended  through  the  outer 

|  steel  plates  of  the  fixture  (Figure  4-8). 

The  procedure  used  to  mount  test  coupons  in  the  mounting  fixture  was  as  follows: 

(1)  Place  inner  (knurled)  mounting  blocks  in  the  recess  of  the  aluminum  setup 
template  (Figure  4-9). 

I  (2)  Place  specimen  between  the  inner  (knurled)  mounting  blocks  and  torque  bolts 

to  30  ft-lbs  (Figure  4-9). 


(3)  Place  the  specimen  and  inner  mount  blocks  within  the  outer  steel  plates  of  the 
mounting  fixture  and  secure  the  outer  plates  with  bolts  (Figure  4-8). 


Once  the  fixture  had  been  assembled,  it  was  installed  in  the  Instron  by  removing  the 
grips  (see  Figure  4-6)  and  securing  the  fixture  by  two  steel  pins.  As  one  might  expect,  the 
assembled  fixture  was  heavy  and  cumbersome.  Therefore,  extreme  care  had  to  be  taken  to 
insure  the  specimens  were  not  damaged  when  mounting  the  fixture  in  the  Instron. 

Due  to  the  very  time  consuming  procedure  required  for  testing  specimens  with  the 
mounting  fixture,  the  only  values  of  eccentricity  tested  with  the  mounting  fixture  were 
eccentricities  of  0.3  inches. 

E.  Photoelastic  Testing. 

In  order  to  verify  the  contour  plotting,  which  was  used  to  study  the  effects  of 
eccentricity  (see  Section  in  E),  photoelastic  testing  was  conducted.  Therefore,  photo¬ 
elasticity  was  utilized  to  verify  "gross"  stress  states  predicted  by  the  analytical  study.  For 
completeness,  a  brief  theory  of  photoelasticity  is  presented  in  Appendix  E. 

Since  the  analytical  study  (Section  IH  E)  indicated  that  the  use  of  the  mounting  fixture 
would  not  appreciably  affect  the  stress  state  adjacent  to  the  discontinuity,  the  mounting 
fixture  was  not  used  during  photoelastic  testing.  Furthermore,  since  a  limited  number  of 
Gr/PEEK  panels  were  available  for  testing,  only  four  specimens  were  tested  using 
photoelasticity.  These  specimens  had  an  eccentricity  of  zero  (see  Figure  4-2),  and  one 
specimen  from  each  of  the  four  ply  lay-ups  was  tested.  The  zero  eccentricity  specimens 
were  used  because  these  specimens  were  available  as  project  residue  from  testing 
completed  by  [10]. 

A  photoelastic  coating  was  bonded  to  each  test  specimen  as  shown  in  Figure  4- 10.  The 
specimens  were  then  mounted  in  the  Instron  and  loaded  at  a  constant  crosshead  velocity  of 
0.05  inches  per  minute.  This  loading  was  continued  until  the  specimen  failed. 

The  photoelastic  testing  was  recorded  on  video  tape,  and  photographs  were  later 
obtained  from  this  tape  (see  Section  V  A).  The  voltmeter  was  also  used  to  obtain  the 


4-12 


magnitude  of  the  applied  loading  at  any  time.  The  voltmeter  was  an  essential  part  of  the  test 
setup  since  a  visible  indication  of  the  applied  load  was  required  to  correlate  the  level  of 
loading  with  the  photographs  obtained  from  the  video  tape  (see  Section  V  A).  The  volt¬ 
meter  is  shown  to  the  left  of  the  specimen  (Figure  4-11). 


.  .  indicates  photoelastic  coating  (Photolastic,  Inc.  PS  1C  coating, 

.044"  thick,  PCI  glue  used  to  bond  the  coating  to  the  specimen) 


Figure  4-10.  Application  of  Photoelastic  Coating 


In  addition  to  the  Instron,  video  camera,  and  voltmeter,  a  polariscope  [  16]  was  used 
during  testing.  Simply  stated,  the  polariscope  (for  purposes  of  this  study)  was  used  to 
polarize  the  light  source  used  during  testing.  Figure  4-11  shows  a  view  (looking  through 
the  video  camera)  of  the  test  setup.  This  figure  shows  the  light  source  on  the  far  right  side 
of  the  photograph.  This  light  source  was  passed  through  a  polarizer  and  then  entered  the 
photoelastic  coating  (center  of  Figure  4-11).  The  light  then  reflected  off  the  reflective 
backing  of  the  photoelastic  coating  (see  Figure  4-10),  and  passed  through  another  polarizer 
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(left  side  of  Figure  4-11).  The  visible  fringes  which  resulted  from  this  process  were  then 
recorded  by  the  video  camera. 


Figure  4-11.  Photoelastic  Test  Setup 

A  simple  schematic  is  shown  in  Figure  4-12.  This  figure  shows  pictorially  the 
photoelastic  technique  used  in  this  study. 

F.  Test  Summary. 

A  total  of  fifty-two  specimens  were  tested  experimentally.  A  test  summary  is  shown  in 
Tables  4-1  through  4-4.  These  tables  are  grouped  according  to  the  ply  lay-up  of  the 
laminates,  and  specify  whether  the  particular  test  specimen  was  instrumented  with  strain 
gages  or  was  used  for  photoelastic  testing.  Finally,  the  tables  show  the  value  of  eccentri¬ 
city  of  the  particular  specimen  and  whether  or  not  the  mounting  fixture  was  utilized. 
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Figure  4-12.  Photoelastic  Technique 


Coupon  Number  Serial  Number  Eccentricity  Strain  Gages  Photoelasticity  Fixture 


Coupon  Number  Serial  Number  Eccentricity  Strain  Gages  Photoelasticity _ Fixture 

B5-90T-00-01-H  014  0.0"  V 

B6-90T-00-02H  015  (HT  1 

B6-90T-00-03-H  016  O0"  V 


Table  4-2.  Test  Summary  for  [90.,  ]  Laminates. 


Table  4-4.  Test  Summary  for  [0° /±45°  I9&  ]2s  Laminates. 


V.  Results  and  Discussion 


This  section  discusses  the  results  from  both  the  analytical  (Section  HI)  and  the  experi¬ 
mental  (Section  IV)  portions  of  this  study.  As  previously  mentioned  in  Section  IV,  ex¬ 
perimentation  was  conducted  for  two  reasons.  First,  to  verify  the  validity  of  the  analytic 
stress  state  predictions  of  the  linear  finite  element  program.  And  secondly,  to  verify  the 
nonlinear  finite  element  predictions  of  ultimate  strength,  failure  modes,  and  stress-strain 
response  of  the  Gr/PEEK  laminates  studied. 

Specifically,  this  section  discusses  the  results  of  both  the  analytical  and  experimental 
studies  to  include: 

(A)  Stress  States 

(B)  Failure  loads/Failure  Modes/Stress-Strain  Response 

(C)  Eccentricity/Mounting  Fixture  -  -  Effects 

(A)  Stress  States. 

The  validity  of  the  analytic  stress  state  predictions  obtained  from  the  linear  finite  element 
program,  were  verified  through  the  use  of  photoelastic  techniques  (Section  TV  E).  Contour 
plots  were  generated  from  data  obtained  from  the  linear  finite  element  program  (Section  HI 
E  (1)).  These  contour  plots,  as  well  as  photographs  taken  during  photoelastic  testing,  are 
shown  in  Figures  5-1  through  5-4.  Each  of  these  figures  corresponds  to  a  particular 
stacking  sequence  and  also  to  a  specific  level  of  loading.  As  already  noted  in  Section  IV  E, 
photoelasticity  was  used  only  to  verify  "gross"  stress  state  predictions  of  the  linear  finite 
element  program. 

Close  inspection  of  Figures  5-1  through  5-4  reveal  that  the  photoelastic  fringes  are  not 
geometrically  symmetric.  Since  the  specimens  tested  using  photoelasticity  were  zero 
eccentricity  specimens  (see  Figure  4-2),  nonsymmetry  of  the  fringes  was  not  expected. 
However,  this  nonsymmetry  can  be  attributed  to  several  different  factors.  First,  photo¬ 
elasticity  (adjacent  to  a  discontinuity)  does  not  provide  an  accurate  measure  of  the  state  of 
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applied  loading=5508  lbs  /  contour  plot 


applied  loading=5508  lbs  /photoelastic  fringes 


Figure  5-1.  Analytical  versus  Experimental  Contours  ( [0,6]  laminate) 
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applied  loading=367.2  lbs  /  contour  plot 


stress  [16].  This  is  due  primarily  to  edge  effects  present  where  the  photoelastic  coating 
meets  the  edge  of  the  discontinuity.  Secondly,  it  is  impossible  to  obtain  a  "perfect" 
adhesive  bond  (i.e.  no  voids,  uniform  glue  thickness,  etc.)  between  the  Gr/PEEK 
specimens  and  the  photoelastic  coatings  (see  Figure  4-10).  Therefore,  the  stress  state  of 
the  Gr/PEEK  specimens  could  not  be  exactly  duplicated  by  the  photoelastic  coating  (see 
Appendix  E  for  a  brief  theory  of  photoelasticity).  Finally,  it  is  impossible  to  locate  a 
circular  discontinuity  in  the  exact  geometric  center  of  a  specimen.  Even  the  slightest  offset 
of  the  discontinuity  will  cause  nonsymmetry  of  the  fringes. 

Figures  5-1  through  5-4  show  good  conelation  between  the  analytical  predictions  of  the 
"gross"  stress  states  and  the  actual  stress  states  present  during  experimentation.  This 
correlation  (for  all  four  stacking  sequences)  is  significant.  As  previously  discussed  in 
Section  II  A(2),  an  analytical  attempt  was  made  to  study  the  Gr/PEEK  laminates  with  a 
two-dimensional  model.  With  this  assumption,  one  assumes  that  through  the  thickness 
normal  stresses  (as  well  as  other  interlaminar  stresses)  may  be  neglected.  Therefore,  since 
the  analytical  predictions  of  the  "gross"  stress  states  correlated  well  with  the  experimental 
stress  states,  it  appears  that  the  two-dimensional  model  adequately  predicts  the  "gross" 
stress  states  observed  from  the  photoelastic  testing.  This  was  important  since  the  linear 
finite  element  program  was  used  to  develop  the  finite  element  model  used  for  this  study 
(see  Section  HI). 

(B)  Failure  Loads/Failure  Modes/Stress-Strain  Response. 

In  this  section,  the  experimental  and  analytical  (nonlinear)  predictions  of  the  failure 
characteristics  of  the  Gr/PEEK  laminates  are  discussed.  First,  the  experimental  failure 
loads,  experimental  scatter,  and  the  analytical  predictions  of  failure  loads  are  presented. 
Secondly,  the  eccentricity  effects  on  the  analytical  predictions  of  the  nominal  (far  field) 
failure  stress  are  discussed.  And  finally,  the  experimental  and  analytical  predictions  of  the 
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failure  modes  and  stress-strain  responses  are  discussed.  This  discussion  is  presented 
according  to  the  stacking  sequence  of  the  laminates  as  follows: 

(1)  [0°  id  laminates 

(2)  [90°  15]  laminates 

(3)  [±45°]4S  laminates 

(4)  t0°/±45o/90o]2s  laminates 

Since  eccentricity  and  boundary  condition  effects  were  observed  to  cause  similar 
analytical  results  for  all  four  of  the  stacking  sequences,  Section  V  C  presents  a  general 
discussion  on  eccentricity  and  boundary  condition  effects  for  all  four  stacking  sequences. 

(i)  [0°16]  laminates- 

Table  5-1  shows  the  experimental  results  of  the  ultimate  strength  tests  for  the  [O0^ 
laminates.  The  scatter  of  the  experimental  data  (Table  5-1)  is  shown  graphically  in  Figure 
5-5.  This  figure  also  lists  the  standard  deviation  for  each  hole  location  (value  of 
eccentricity)  tested  experimentally.  The  standard  deviation  was  calculated  as: 


1 

2 


r  n 


TT 


(l) 


where: 

s  =  standard  deviation 
y,  =  i*  data  point 
N  =  total  number  of  data  points 
y  =  mean  value  of  all  data  points 
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Serial  Number 

Eccentricity 

Mounting 

Fixture 

Failure  Load 
(lbs) 

002 

0.0" 

4918.12 

003 

0.0" 

5174.97 

004 

0.0" 

5376.14 

005 

0.1" 

6987.95 

006 

0.1" 

6675.78 

007 

0.1" 

7090.03 

008 

0.3" 

9545.49 

009 

0.3" 

9421.02 

010 

0.3" 

9468.81 

Oil 

0.3" 

✓* 

6986.00 

012 

0.3" 

V 

9085.47 

013 

0.3" 

V 

8345.11 

*  (Vindicates  that  the  mounting  fixture  was  used) 

Figure  5-5  shows  that  for  the  eccentricity=0.3"  testing  (with  the  mounting  fixture), 
there  was  an  unacceptable  amount  of  experimental  scatter.  This  scatter  may  have  been  a 
result  of  initial  damage  to  the  specimens  (note  the  relatively  low  failure  load  for  specimen 
Oil,  Table  5-1).  As  noted  in  Section  IV  D,  testing  with  the  mounting  fixture  was  difficult 
due  to  the  weight  of  the  assembled  fixture.  Therefore,  it  would  have  been  very  easy  to 
damage  a  specimen  while  installing  the  mounting  fixture  in  the  Instron  Test  Machine. 

A  summary  of  the  average  experimental  failure  loads  (calculated  from  Table  5-1),  as 
well  as  the  analytical  predictions  of  failure  loads,  are  given  in  Table  5-2.  In  addition.  Table 
5-2  shows  the  percent  difference  between  the  average  experimental  failure  loads  and  the 
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Figure  5-5.  Experimental  Scatter  for  [0°6  ]  laminates 


analytical  predictions  of  failure  loads.  This  percent  difference  was  calculated  as: 


avg  experimental  failure  load  (lbs)  -  analytical  failure  load  (lbs) 
avg  experimental  failure  load  (lbs) 


x  100 


If  the  analytical  predictions  of  the  nominal  (far  field)  failure  stresses  are  examined,  the 
effects  of  eccentricity  on  the  nominal  failure  stress  can  be  determined.  Table  5-3  lists  both 
the  analytical  failure  loads  and  the  nominal  (far  field)  failure  stresses. 
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Table  5-2.  Summary  of  Average  Experimental  Failure  Loads  and 
Analytical  Failure  Loads  for  [0°6  ]  Laminates. 

Eccentricity 

Average  Experimental 
Failure  Load  (lbs) 

Analytical  Failure 
Load  (lbs) 

%  diff. 

(Ex.  vs  Anal.) 

0.0" 

(no  fixture) 

5156.41 

5935.36 

-15.1% 

0.1" 

(no  fixture) 

6917.92 

7041.71 

-1.8% 

0.3" 

(no  fixture) 

9478.44 

9241.27 

+2.5% 

0.3" 

(with  fixture) 

8138.86 

8952.14 

-9.1% 

Table  5-3.  Summary  of  Analytical  Failure  Loads  and 

Nominal  (far  field)  Failure  Stresses  for  [0°16  ]  Laminates. 

Eccentricity 

Analytical  Failure 

Load  (lbs) 

Nominal  Failure  (far 
field)  Stress  (psi) 

0.0" 

(no  fixture) 

5935.36 

58,882.54 

0.1" 

(no  fixture) 

7041.71 

59,878.49 

0.3" 

(no  fixture) 

9241.27 

61,119.51 

0.3" 

(with  fixture) 

8952.14 

59,207.27 

The  initial  attempt  in  this  thesis  was  to  remove  the  edge  effects  created  by  locating  the 
discontinuity  adjacent  to  the  edge  of  the  specimen.  Therefore,  for  all  the  specimens 
considered  in  this  study,  the  distance  from  the  edge  of  the  discontinuity,  to  the  edge  of  the 
specimen  (see  Figure  4-2),  was  always  greater  than  or  equal  to  the  diameter  of  the  discon- 
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tinuity  (0.4").  This  distance  was  maintained  (for  all  specimens)  so  that  edge  effects  would 
not  vary  for  the  different  values  of  eccentricity  considered  in  this  study.  From  Table  5-3, 
note  that  the  analytical  predictions  of  nominal  failure  stresses  do  not  appreciably  differ  for 
the  different  values  of  eccentricity;  therefore,  since  edge  effects  were  effectively  controlled 
by  main-taining  at  least  0.4  inches  between  the  edge  of  the  discontinuity  and  the  edge  of 
the  specimen,  it  appears  that  the  values  of  eccentricity  considered  in  this  study  did  not 
appreciably  affect  the  nominal  failure  stress  for  the  [O^g]  laminates. 

The  "knock  down"  (reduction)  of  the  nominal  failure  stress,  due  to  the  discontinuity, 
was  significant.  The  nominal  failure  stress  for  [0°ig]  laminates  (without  a  discontinuity)  is 
298,300  psi  [13].  Therefore,  for  this  laminate,  the  nominal  failure  stress  (due  to  the  0.4" 
discontinuity)  was  reduced  by  a  factor  of  approximately  5  (from  Table  5-3).  This  is 
significandy  higher  than  the  knockdown  for  an  isotropic  material  of  the  same  geometry. 
Peterson  [15]  reports  that  for  an  isotropic  material  (of  the  same  geometry),  the  nominal 
failure  stress  would  be  reduced  by  a  factor  of  3.46.  Since  the  [0°16]  laminates  failed  along 
the  fiber  direction,  the  magnitude  of  the  knockdown  of  nominal  failure  stress  is  not  un¬ 
reasonable.  This  is  because  once  the  ultimate  stress  of  the  laminate  was  reached  (at  the  edge 
of  the  discontinuity),  the  laminate  split  along  the  fiber  direction.  This  failure  mode  occurred 
at  a  relatively  low  nominal  stress.  Once  the  zero  degree  fibers  adjacent  to  the  discontinuity 
failed,  the  matrix  (which  is  relatively  weak)  failed  very  rapidly  (split). 

For  all  of  the  [0°  15]  specimens,  the  failure  mode  was  characterized  by  splitting  at  the 
edge  of  the  discontinuity  (parallel  to  the  fibers).  Figures  5-6  and  5-7  show  photographs  of 
the  failed  specimens.  The  splitting  of  the  specimens  (cracks)  has  been  highlighted  by  silver 
lines  (Figures  5-6  and  5-7).  The  terms  "hinge"  and  "no  hinge"  shown  on  these  figures, 
refers  to  testing  with  or  without  the  mounting  fixture,  respectively.  For  purposes  of  this 
study,  failure  (both  experimental  and  analytical)  was  assumed  to  occur  at  the  load  which 
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produced  splitting.  This  assumption  was  made  because  once  the  specimens  split,  the 
geometry  of  the  problem  was  entirely  changed.  Instead  of  a  specimen  with  a  circular 
discontinuity,  the  specimen  was  now  comprised  of  three  different  sections  (see  Figure  5-6). 
These  three  sections  could  be  described  as  two  rectangular  sections  above  and  below  the 
discontinuity,  and  a  center  section  (left  and  right  of  the  discontinuity)  which  was  unloaded. 
Therefore,  the  original  finite  element  models  (Figure  3-3)  could  not  be  expected  to  provide 
accurate  data  beyond  the  "splitting"  load.  This  is  due  to  the  fact  that  the  geometry  of  the 
finite  element  model  no  longer  matched  the  geometry  of  the  "split"  specimens. 

The  determination  of  the  experimental  failure  load  was  an  easy  task.  When  splitting 
occurred,  an  instantaneous  drop  in  the  applied  load  was  observed  on  the  Instron  strip  chart 
(see  Figure  4-6).  Also,  the  output  strain  gage  data  (Section  IV  C)  showed  discontinuous 
jumps  in  strain  at  the  instant  of  splitting.  Therefore,  an  accurate  value  of  the  experimental 
failure  load  was  easily  obtained. 

Once  the  experimental  failure  mode  had  been  observed,  the  determination  of  the 
analytical  failure  load  was  also  straightforward.  Since  splitting  was  initiated  at  the  edge  of 
the  discontinuity  (point  A,  Figure  5-8),  the  stress  at  the  edge  of  the  discontinuity  was 
obtained  (from  the  nonlinear  output)  for  each  increment  of  loading.  Since  the  nonlinear 
finite  element  program  provided  stress  and  strain  data  only  at  the  centroid  of  the  finite 
elements,  Lagrange  extrapolation  [5]  was  used  to  calculate  the  value  of  stress  at  the  edge  of 
the  discontinuity.  Once  the  stress  level  at  the  edge  of  the  discontinuity  reached  the  ultimate 
stress  of  the  Gr/PEEK,  splitting  was  assumed  to  occur,  and  the  failure  load  from  the 
nonlinear  output  was  recorded. 

One  interesting  point  regarding  the  splitting  of  these  laminates  becomes  apparent  if  the 
values  of  transverse  stress  (oy)  along  the  split  (points  B,C,and  D;  Figure  5-8)  are  examined. 

It  was  observed  that  these  transverse  stresses  (less  than  200  psi)  were  too  small  to  cause 
splitting.  Therefore,  the  splitting  could  be  caused  by  a  dynamic  process,  such  as  a  relax- 
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Figure  5-8.  [0°16  ]  Laminate  Splitting 


ation  stress  wave  generated  by  splitting  of  the  fibers  at  the  hole.  Therefore,  if  splitting  was 
caused  by  a  dynamic  phenomenon  (such  as  a  stress  wave),  the  finite  element  solution  (a 
static  solution)  can  not  be  expected  to  predict  splitting.  Other  possibilities  for  splitting  might 
include  edge  effects  or  some  other  three-dimensional  effect  at  the  edge  of  the  discontinuity. 

Stress-strain  response  plots  were  constructed  by  [13]  for  symmetric  models  (eccentri¬ 
city  =  0.0").  Therefore,  in  an  effort  to  bound  the  ranges  of  eccentricities  considered  in  this 
study,  only  the  stress-strain  response  plots  for  the  eccentricity  =  0.3"  models  were  con¬ 
structed  (Figures  5-9  through  5-1 1).  Figure  5-9  shows  the  experimental  and  analytical  near 
field  response  at  strain  gage  #1.  The  experimental  and  analytical  near  field  response  at  strain 
gage  #2  is  shown  in  Figure  5-10.  Finally,  the  experimental  and  analytical  far  field  response 
is  shown  in  Figure  5-1 1  as  the  average  strains  from  strain  gages  #3  and  #4.  From  these 
three  figures,  note  that  there  is  excellent  (within  5%)  correlation  between  the  analytical  pred¬ 
ictions  of  axial  stress-strain  response,  and  the  experimental  axial  stress-strain  response  (for 
a  constant  value  of  stress).  This  correlation  was  noted  at  both  the  near  field  (Figures  5-9 
and  5-10)  and  the  far  field  (Figure  5-11)  locations.  Also,  from  these  three  figures,  note  that 
there  is  good  (within  10%)  correlation  between  the  analytical  predictions  of  transverse 
stress-strain  response  (for  levels  of  stress  up  to  the  point  of  the  splitting  of  the  laminate)  and 
the  experimental  transverse  stress-strain  response  (for  a  constant  value  of  stress). 
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Figure  5-9.  Near  Field  (gage  #1)  Stress-Strain  Response 
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1.  *  indicates  serial  number  of  test  specimen 

2.  axial  strain  for  specimen  #010  not  available  (bad  gage) 

3.  0  degree  lay-up  /  eccentricity=0.3"  /  no  mounting  fixture 


Figure  5-10.  Near  Field  (gage  #2)  Stress  Strain  Response 
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100000 


STRAIN  GAGE  LOCATIONS 
gage#4  (reverse  side) 

m 

gag<#2-g  ym 
gage  #1— #  ,P##3 

H 

NOTES: 

1.  *  indicates  serial  number  of  test  specimen 

2.  0  degree  lay-up  /  eccentricity=0.3"  /  no  mounting  fixture 


Figure  5-11.  Far  Field  (gage  #3)  Stress  Strain  Response 


In  summary,  the  analytical  solution  did  a  good  job  of  predicting  the  splitting  (failure) 
load  of  the  [0°16]  laminates.  Also,  the  analytical  stress-strain  response  predictions 
correlated  well  with  the  experimental  stress-strain  response.  Finally,  the  analytical  solution 
did  not  predict  the  failure  progression  (splitting)  of  these  laminates.  As  already  mentioned, 
this  splitting  may  be  caused  by  a  dynamic  phenomenon  (such  as  a  stress  wave),  or  some 
other  three-dimensional  phenomenon  (such  as  edge  effects).  Therefore,  the  two-dimen¬ 
sional  finite  element  solution  (a  static  solution)  was  not  able  to  predict  splitting. 

(2)  [90°16]  laminates- 

Table  5-4  shows  the  experimental  results  of  the  ultimate  strength  tests  for  the  [90°  i&] 
laminates.  The  scatter  of  the  experimental  data  (Table  5-4)  is  shown  graphically  in  Figure 
5-12.  This  figure  also  lists  the  standard  deviation  for  each  hole  location  (value  of  eccen¬ 
tricity)  tested  experimentally.  The  standard  deviation  was  calculated  from  equation  (1), 
Section  V  B(l).  Figure  5-12  shows  that  experimental  scatter  was  not  significant  during  the 
testing  of  these  laminates. 

A  summary  of  the  average  experimental  failure  loads  (calculated  from  Table  5-4)  as  well 
as  the  analytical  predictions  of  failure  loads  are  given  in  Table  5-5.  In  addition.  Table  5-5 
shows  the  percent  difference  between  the  average  experimental  failure  loads  and  the  analy¬ 
tical  predictions  of  failure  loads.  This  percent  difference  was  calculated  from  equation  2, 
Section  V  B(l). 

If  the  analytical  predictions  of  the  nominal  (far  field)  failure  stresses  are  examined,  the 
effects  of  eccentricity  on  the  nominal  failure  stress  can  be  determined.  Table  5-6  lists  both 
the  analytical  failure  loads  and  the  nominal  (far  field)  failure  stresses.  As  previously 
discussed  (Section  V  B(l)),  the  distance  from  the  edge  of  the  discontinuity  to  the  edge  of 
the  specimen  (Figure  4-2)  was  always  greater  than  or  equal  to  the  diameter  of  the  discon¬ 
tinuity  (0.4").  This  distance  was  maintained  so  that  edge  effects  would  not  vary  for  the 
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Table  5-4.  Experimental  Results  for  [90  jg  ]  Laminates. 


Serial  Number 

Eccentricity 

015 

0.0" 

016 

0.0" 

017 

0.0" 

018 

0.1" 

019 

0.1" 

020 

0.1" 

021 

0.3" 

022 

0.3" 

023 

0.3" 

024 

0.3" 

025 

0.3" 

026 

0.3" 

Mounting 

Fixture 


Failure  Load 
(lbs) 


704.40 


742.29 


799.78 


845.85 


828.64 


956.03 


1085.12 


1198.66 


1247.73 


1213.36 


1301.00 


1150.74 


*<✓  indicates  that  the  mounting  fixture  was  used) 


different  values  of  eccentricity  considered  in  this  study.  From  Table  5-6,  note  that  the 
analytical  predictions  of  nominal  failure  stresses  do  not  appreciably  differ  for  the  different 
values  of  eccentricity;  therefore,  since  edge  effects  were  controlled  by  maintaining  at  least 
0.4  inches  between  the  edge  of  the  discontinuity  and  the  edge  of  the  specimen,  it  appears 
that  the  values  of  eccentricity  considered  in  this  study  did  not  appreciably  affect  the  nominal 
failure  stress  for  this  stacking  sequence. 

The  knock  down  (reduction)  of  the  nominal  failure  stress  (due  to  the  0.4"  discontinuity) 
was  significantly  lower  than  that  obtained  for  the  [0°ig]  laminates.  The  nominal  failure 
stress  for  [90°  J6]  laminates  (without  a  discontinuity)  is  14,1 10  psi  [13].  Therefore,  for 
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Figure  5-12.  Experimental  Scatter  for  [90  16  ]  Laminates 


Std  Dev 
61.6  lbs  * 


1  Std  Dev  T 
\  39.2  lbs  i 


Std  Dev=standard  deviation 

*  indicates  that  the  mounting  fixture  was  used 

•  indicates  the  average  failure  load  (mean) 


Std  Dev 
56.4  lbs 


Std  Dev  I 
68.1  lbs  I 


Table  5-5.  Summary  of  Average  Experimental  Failure  Loads  and 
Analytical  Failure  Loads  for[90°6  ]  Laminates. 

Eccentricity 

Average  Experimental 
Failure  Load  fibs) 

Analytical  Failure 
Load  fibs) 

%  diff. 

(Ex.  vs  Anal.) 

0.0" 

(no  fixture) 

748.82 

729.05 

+2.6% 

0.1" 

(no  fixture) 

876.84 

894.89 

m 

0.3" 

(no  fixture) 

1177.17 

1217.37 

-3.4% 

0.3" 

(with  fixture) 

1221.70 

1172.56 

+4.0% 
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Table  5-6.  Summary  of  Analytical  Failure  Loads  and 

Nominal  (far  field)  Failure  Stresses  for  [90  °6  ]  Laminates. 

Eccentricity 

Analytical  Failure 

Load  (lbs) 

Nominal  Failure  (far 
field)  Stress  (psi) 

0.0” 

(no  fixture) 

729.05 

7,232.64 

0.1” 

(no  fixture) 

894.89 

7,609.61 

0.3" 

(no  fixture) 

1217.37 

8,051.39 

0.3" 

(with  fixture) 

1172.56 

7,755.03 

this  laminate,  the  nominal  failure  stress  (due  to  the  0.4"  discontinuity)  was  reduced  by  a 
factor  of  approximately  1.8  (see  nominal  failure  stresses.  Table  5-6).  This  is  significantly 
lower  than  the  knock  down  of  nominal  failure  stress  (due  to  a  circular  discontinuity)  for  an 
isotropic  material  of  the  same  geometry,  reported  as  3.46  by  [15]. 

For  all  of  these  specimens,  experimental  failure  occurred  parallel  to  the  fibers  and  was 
initiated  at  the  discontinuity.  Figures  5-13  and  5-14  show  photographs  of  the  failed 
specimens.  This  mode  of  failure  was  expected  since  the  [90oig]  laminates  are  charact¬ 
erized  by  matrix  failure  (parallel  to  the  fibers)  at  the  narrowest  part  of  the  cross-sectional 
area  [13].  Failure  of  these  laminates  was  very  rapid,  therefore,  it  was  impossible  to 
determine  if  failure  was  initiated  above  or  below  the  discontinuity. 

Figures  5-15  through  5-18  show  the  analytical  predictions  of  failure  progression  and 
failure  loads  for  the  [90°ig]  laminates.  These  figures  show  that  failure  of  these  laminates 
was  initiated  at  the  discontinuity.  These  figures  also  show  that  "free-floating"  nodes 
occurred  at  the  analytical  failure  load.  The  term  "free-floating"  node  indicates  an  instability 
of  the  displacements  of  a  node.  This  instability  occures  due  to  the  large  displacement  of  a 
node,  which  causes  the  rotation  of  a  line  connecting  two  nodes  (i.e.  one  side  of  a  finite 
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I 

I 


Figure  5-13.  Failed  [90^  ] Specimens  (no  mounting  fixture) 


Figure  5-14.  Failed  [90  £5  ]  Specimen  (mounting  fixture) 
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The  shaded  portion  of  the  model  is 
magnified  below,  and  the  displacements 
of  tins  section  are  magnified  100  times. 


E)  indicates  element  failure  at  688.4  lbs 


1.  Free-floating  nodes  occurred  (instability)  @  729.05  lbs;  in  addition,  the 
load  increments  converged  at  the  FAILURE  LOAD  (729.05  lbs) 


Figure  5-15.  Failure  Progression  (90  degree  /  e*G.0"  /  no  mounting  fixture) 
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■  indicates  element  failure  at  8 1 4.73  lbs 

B  indicates  element  failure  at  834.02  lbs 

□  indicates  element  failure  at  865.09  lbs 

■  indicates  element  failure  at  894.89  lbs  (FAILURE  LOAD) 


1.  Free-floating  nodes  occurred  (instability)  @  894.89  lbs;  in  addition,  the 
load  increments  converged  at  the  FAILURE  LOAD  (894.89  lbs). 


Figure  3-16.  Failure  Progression  (90 degree  /e*0.1"  /no  mounting  fixture) 
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The  shaded  portion  of  the  model  is 
magnified  below,  and  the  displacements 
of  this  section  are  magnified  100  times. 


■  indicates  element  failure  at  1062.09  lbs 
E3  indicates  element  failure  at  1 129.33  lbs 
D  indicates  element  failure  at  1 179.75  lbs 
©  indicates  element  failure  at  1217.37  lbs 


1.  Free-floating  nodes  occurred  (instability)  @  1 179.75  lbs;  in  addition,  the 
load  increments  converged  at  the  FAILURE  LOAD  (1217.37  lbs). 


Figure  5- 17.  Failure  Progression  (90  degree  /  e=0.3"  /  no  mounting  fixture) 
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The  shaded  portion  of  the  model  is 


magnified  below,  and  the  displacements 


a  indicates  element  failure  at  1042.59  lbs 

d  indicates  element  failure  at  1 1 10.76  lbs 

O  indicates  element  failure  at  1 120.92  lbs 

•  indicates  element  failure  at  1 172.56  lbs 


1.  Free-floating  nodes  occurred  (instability)  @  1 172.56  lbs;  in  addition,  the 
load  increments  converged  at  the  FAILURE  LOAD  (1 172.56  lbs). 


element).  This  rotation  of  a  line,  connecting  two  nodes,  is  visible  in  Figures  5-15  through 
5-18. 

Stress-strain  response  plots  for  the  eccentricity  =  0.3"  model  were  constructed  (Figures 
5-19  through  5-21).  Figure  5-19  shows  die  experimental  and  analytical  near  field  re¬ 
sponse  at  strain  gage  #1.  The  experimental  and  analytical  near  field  response  at  strain  gage 
#2  is  shown  in  Figure  5-20.  And  finally,  the  experimental  and  analytical  far  field  response 
is  shown  in  Figure  5-21  as  the  average  strains  from  strain  gages  #3  and  #4.  From  these 
three  figures,  note  that  there  is  good,  within  10%,  correlation  between  the  analytical 
predictions  of  axial  stress-strain  response,  and  the  experimental  axial  stress-strain  response 
(for  a  constant  value  of  stress).  Also,  from  these  figures,  note  that  there  is  excellent 
correlation  between  the  analytical  predictions  of  transverse  stress-strain  response  and  the 
experimental  transverse  stress-strain  response  (for  a  constant  value  of  stress). 

In  summary,  the  analytical  solution  did  an  excellent  job  of  predicting  both  the  failure 
loads  and  failure  modes  for  the  [90Oj6]  laminates.  In  addition,  the  analytical  stress-strain  _ 
response  predictions  correlated  very  well  with  the  experimental  stress-strain  response. 

(3)  [±45°]4s  laminates- 

Table  5-7  shows  the  experimental  results  of  the  ultimate  strength  tests  for  the  [±45°]4s 
laminates.  The  scatter  of  the  experimental  data  (Table  5-7)  is  shown  graphically  in  Figure 
5-22.  This  figure  also  lists  the  standard  deviation  for  each  hole  location  (value  of  eccen¬ 
tricity)  tested  experimentally.  The  standard  deviation  was  calculated  from  equation  (1), 
Section  V  B(l).  Figure  5-22  shows  that  experimental  scatter  was  not  significant  during  the 
testing  of  these  laminates. 

A  summary  of  the  average  experimental  failure  loads  (calculated  from  Table  5-7),  as  well 
as  the  analytical  predictions  of  failure  loads,  are  given  in  Table  5-8.  In  addition,  Table  5-8 


Applied  Stress  (psi) 


STRAIN  GAGE  LOCATIONS 
gage#4  (reverse  side) 

Wa 

gage#2— •  ^ 
gage#l— •  gage#3 

NOTES: 

1.  *  indicates  serial  number  of  test  specimen 

2.  90  degree  lay-up  /  eccentricity=0.3"  /  no  mounting  fixture 


Figure  5-19.  Near  Field  (gage  #1)  Stress  Strain  Response 
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Applied  Stress  (psi) 
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pplied  Stress  (psi) 


Strain  (micro  in/in) 


STRAIN  GAGE  LOCATIONS 
gage#4  (reverse  side) 

1 

gage#2~o  /* 
gage#l  — •  S»ge#3 

NOTES: 

1.  *  indicates  serial  number  of  test  specimen 

2.  90  degree  lay-up  /  eccentricity=0.3"  /  no  mounting  fixture 


Figure  5-21.  Far  Field  (gage  #3)  Stress  Strain  Response 
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Table  5-7.  Experimental  Results  for  [±45°]  ^Laminates. 

Serial  Number 

Eccentricity 

Mounting 

Fixture 

Failure  Load 
(lbs) 

028 

0.0" 

2434.54 

029 

0.0" 

2289.37 

030 

0.0" 

2290.32 

031 

0.1" 

2807.02 

032 

0.1" 

3044.05 

033 

0.1" 

3261.84 

034 

0.3" 

4122.58 

035 

0.3" 

4030.70 

036 

0.3" 

3813.52 

037 

0.3" 

4537.12 

038 

0.3" 

v/ 

4421.36 

039 

0.3" 

✓ 

4105.35 

*  ( indicates  that  the  mounting  fixture  was  used) 

shows  the  percent  difference  between  the  average  experimental  failure  loads  and  the  analy¬ 
tical  predictions  of  failure  loads.  This  percent  difference  was  calculated  from  equation  (2), 
Section  V  B(l). 

The  failure  loads  predicted  analytically  were  significantly  lower  than  the  experimental 
failure  loads  (see  Table  5-8).  This  difference  may  be  caused  by  a  phenomenon  known  as 
scissoring.  This  phenomenon  can  best  be  described  as  a  tendency  of  the  ±45  degree  fibers 
to  align  themselves  with  the  load.  This  phenomenon  was  noted  by  [10,13]  and  is  shown  in 
Figure  5-23.  The  result  of  scissoring  is  that  the  laminate  is  able  to  sustain  higher  loads. 
This  is  due  to  the  fact  that  scissoring  causes  the  fibers,  the  principle  load  carriers,  to 
become  more  aligned  with  the  load  (see  Figure  5-23).  Therefore,  the  laminate  is  able  to 
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StdDev  _ 
129.6  lbs  I 


I  StdDev 
I  182.5  lbs* 


StdDev  I 
185.7  lbs  1 


3  2000 

'3 

U< 

1500 


StdDev 
68.2  lbs 


Std  Dev=standard  deviation 

*  indicates  that  the  mounting  fixture  was  used 

•  indicates  the  average  failure  load  (mean) 


0.0  o.l  0.3 

Eccentricity  (in) 


Figure  5-22.  Experimental  Scatter  for  [±45  ]  .  Laminates 


Table  5-8.  Summary  of  Average  Experimental  Failure  Loads  and 
Analytical  Failure  Loads  for  [±45  °]4sLaminates. 


Eccentricity 

Average  Experimental 
Failure  Load  (lbs) 

Analytical  Failure 
Load  (lbs) 

%  diff. 

(Ex.  vs  Anal.) 

0.0" 

(no  fixture) 

2338.08 

1949.26 

+16.6% 

0.1" 

(no  fixture) 

3037.64 

2133.70 

+29.8% 

0.3" 

(no  fixture) 

3988.93 

2931.39 

+26.5% 

0.3" 

(with  fixture) 

4354.61 

2590.33 

+40.5% 

±  45  degreefibers  45  degrees 


K"  o 


no  applied  load 


Angle  "A" 


LOAD 


-Xo 


load  applied 

Note:  Angle  "A"  is  less  than  45  degrees.  This  decrease  (from  the 
original  45  degrees)  is  caused  by  scissoring.  Also  note  that 
the  fibers  are  tending  to  align  themselves  with  the  loading. 


Figure  5-23.  Scissoring  of  ±45  Degree  Fibers 


sustain  higher  loading.  Since  the  ply  orientations  (±45  degree  is  this  case)  are  fixed  as 
inputs  to  the  nonlinear  material  finite  element  program,  the  nonlinear  analysis  did  not 
account  for  scissoring.  It  should  be  pointed  out  that  scissoring  is  a  nonlinear  strain- 
displacement  phenomenon;  therefore,  it  is  a  geometric  nonlinearity  as  opposed  to  a  material 
nonlinearity  (such  as  stress-strain).  As  a  result,  the  nonlinear  predictions  of  failure  loads 
are  lower  than  those  obtained  experimentally  (Table  5-8).  Another  possibility  for  the 
difference  in  experimental  and  analytical  failure  loads,  is  the  fact  that  the  two-dimensional 
approach  used  in  this  study  did  not  take  into  account  any  of  the  interlaminar  stresses,  or 
other  three-dimensional  effects,  that  may  alter  the  failure  loads  for  this  stacking  sequence. 

If  the  analytical  predictions  of  the  nominal  (far  field)  failure  stresses  are  examined,  the 
effects  of  eccentricity  on  the  nominal  failure  stress  can  be  determined.  Table  5-9  lists  both 
the  analytical  failure  loads  and  the  nominal  (far  field)  failure  stresses.  As  previously  dis¬ 
cussed  (Section  V  B(l)),  the  distance  from  the  edge  of  the  discontinuity  to  the  edge  of  the 
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Tabic  5-9.  Summary  of  Analytical  Failure  Loads  and 
Nominal  (far  field)  Failure  Stresses  for  [±45°]4s  Laminates. 


Eccentricity 

Analytical  Failure 

Load  fibs) 

Nominal  Failure  (far 
field)  Stress  (psi) 

1949.26 

19,337.90 

■2SBM1 

2133.70 

18,143.71 

0.3" 

(no  fixture) 

2931.39 

19,387.50 

0.3" 

(with  fixture) 

2590.33 

17,131.81 

specimen  (Figure  4-2)  was  always  greater  than  or  equal  to  the  diameter  of  the  discontinuity 
(0.4").  This  distance  was  maintained  so  that  edge  effects  would  not  vary  for  the  different 
values  of  eccentricity  considered  in  this  study.  From  Table  5-9,  note  that  the  analytical 
predictions  of  nominal  failure  stresses  do  not  appreciably  differ  for  the  different  values  of 
eccentricity;  therefore,  since  edge  effects  were  controlled  by  maintaining  at  least  0.4  inches 
between  the  edge  of  the  discontinuity  and  the  edge  of  the  specimen,  it  appears  that  the 
values  of  eccentricity  considered  in  this  study  did  not  appreciably  affect  the  nominal  failure 
stress  for  this  stacking  sequence. 

The  knock  down  (reduction)  of  the  nominal  failure  stress,  due  to  the  0.4"  discontinuity, 
was  calculated  from  Table  5-9.  The  nominal  failure  stress  for  these  laminates,  without  a 
discontinuity,  is  equal  to  56,170  psi  [13].  Therefore,  the  nominal  failure  stress  was  re¬ 
duced  by  a  factor  of  approximately  3.  The  reduction  of  nominal  failure  stress,  due  to  a  cir¬ 
cular  discontinuity,  for  an  isotropic  material  of  the  same  geometry,  is  reported  as  3.46  [15]. 

The  failure  mode  of  these  specimens  was  visible  as  a  shear  dominated  failure  surface 
(i.e.  mode  II  dominated  the  failure).  This  shear  dominated  failure  surface  and  the  high 
strains  (visible  as  necking)  near  the  discontinuity,  are  shown  in  Figures  5-24  and  5-25. 


Figure  5-24.  Failed  [±45°]  ^Specimens  (no  mounting  fixture) 


Figure  5-25.  Failed  [±45  °]  ^Specimen  (mounting  fixture) 
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Figures  5-26  through  5*29  show  the  analytical  predictions  of  failure  progression  and 
failure  loads  for  the  [±45°]^  laminates.  From  these  figures,  note  that  the  damage  zones 
(failed  Finite  elements)  appear  to  be  failing  in  ±45  degree  directions.  Therefore,  it  appears 
that  the  analytical  solution  is  predicting  a  shear  (mode  II)  type  failure.  However,  previous 
research  by  [13]  examined  Gr/PEEK  laminates,  for  zero  eccentricity,  which  had  been  tested 
to  a  percentage  of  their  failure  load.  Through  the  use  of  stereo  x-rays,  the  failure  zones 
appeared  as  triangular  areas  above  and  below  the  discontinuity  [13].  These  damaged 
zones,  at  95%  of  failure  load  [13],  are  reproduced  in  Figure  5-30.  The  difference  in 
failure  zones  between  those  predicted  analytically  (this  study)  and  those  shown  in  Figure 
5-30  [13],  might  be  due  to  the  fact  that  interlaminar  stresses,  and  other  three-dimensional 
effects,  were  not  considered  in  this  study. 

Stress-strain  response  plots  for  die  eccentricity  =  0.3"  model  were  constructed  (Figures 
5-31  through  5-33).  Figure  5-31  shows  the  experimental  and  analytical  near  field  response 
at  strain  gage  #1.  The  experimental  and  analytical  near  field  response  at  strain  gage  #2  is 
shown  in  Figure  5-32.  And  finally,  the  experimental  and  analytical  far  field  response  is 
shown  in  Figure  5-33  as  the  average  strains  from  strain  gages  #3  and  #4.  From  these  three 
figures,  note  that  very  little  experimental  data  was  available.  This  is  due  to  the  fact  that  two 
of  the  three  specimens  tested  experimentally  did  not  provide  strain  information.  The  strain 
gages  on  these  two  specimens  "saturated"  at  a  very  low  level  of  load.  This  saturation 
resulted  from  improperly  adjusted  instrumentation  equipment,  and  was  not  the  result  of 
strains  above  the  operating  limits  of  the  gages  (gages  used  were  good  up  to  5%  strain). 

The  strain  gages  on  the  specimen,  which  did  provide  strain  information,  saturated  at 
approximately  4%  far  field  strain  (Figure  5-33),  and  between  1%  and  2%  near  field  strain 
(Figures  5-31  and  5-32).  The  ±45  degree  laminates  are  difficult  (experimentally)  to  obtain 
strain  data  for,  this  is  due  to  high  strains  present  at  the  failure  loads  [10.13].  In  addition, 


The  shaded  portion  of  the  model  is 
magnified  below,  and  all  displacements 
are  magnified  10  times. 


O  indicates  ±45  degree  ply  failures  @  1435.71  lbs 
•  indicates  ±45  degree  ply  failures  @  1448.85  lbs 

□  indicates  ±45  degree  ply  failures  @  1729.56  lbs 
CD  indicates  ±45  degree  ply  failures  @  1757.47  lbs 

□  indicates  ±45  degree  ply  failures  @  1949.26  lbs 

ANALYTICAL  FAILURE  LOAD  WAS  1949.26  LBS 


Figure  5-26.  Failure  Progression  (±45  degree  /  e=0.0"  /  no  mounting  fixture) 
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The  shaded  portion  of  the  model  is 
magnified  below,  and  all  displacements 
are  magnified  10  times. 


■  indicates  ±45  degree  ply  failures  @  1710.67  lbs 
a  indicates  ±45  degree  ply  failures  @  1736. 16  lbs 
q  indicated ±45  degree  ply  failures  @  1920.35  lbs 
•  indicates  ±45  degree  ply  failures  @  1927.32  lbs 
O  indicates  ±45  degree  ply  failures  @  21 17.58  lbs 
O  indicates  ±45  degree  ply  failures  @  2 1 32. 19  lbs 

ANALYTICAL  FAILURE  LOAD  WAS  2132.19  LBS. 


The  shaded  portion  of  the  model  is 
magnified  below,  and  all  displacements 
are  magnified  10  times. 


■  indicates  ±45  degree  ply  failures  @  2403. 13  lbs 
°  indicates  ±45  degree  ply  failures  @  2438.23  lbs 

•  indicates  ±45  degree  ply  failures  @  2652.16  lbs 

•  indicates  ±45  degree  ply  failures  @  2672.68  lbs 
Q  indicates  ±45  degree  ply  failures  @  2733. 1 6  lbs 
H  indicates  ±45  degree  ply  failures  @  2929.94  lbs 


ANALYTICAL  FAILURE  LOAD  WAS  2929.94  LBS. 


Figure  5*28.  Failure  Progression  (±45  degree  /  e=0.3"  /  no  mounting  fixture) 
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The  shaded  portion  of  the  model  is 
magnified  below,  and  all  displacements 
are  magnified  10  times. 


■  indicates  ±45  degree  ply  failures  @  2287.55  lbs 
°  indicates  ±45  degree  ply  failures  @  232 1 .40  lbs 
©  indicates  ±45  degree  ply  failures  @  2549.89  lbs 
•  indicates  ±45  degree  ply  failures  @  2587.96  lbs 

ANALYTICAL  FAILURE  LOAD  WAS  2587.96  LBS. 


Figure  5-29.  Failure  Progression  (±45  degree  /  e=0.3"  /  mounting  fixture) 
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NOTES:  Specimen 

1.  Figure  shows  center  section  of  laminate  only. 

2.  This  drawing  was  created  from  Figure  77  [13]. 

3.  Intent  is  ONLY  to  show  experimentally  obtained  damaged 

areas. 

Figure  5-30.  Damaged  Zones  [±45°]  4s  Laminate  [13] 

scissoring  has  also  been  found  to  contribute  to  the  premature  failure  of  strain  gages  on 
these  laminates  [10,13].  The  analytical  stress-strain  response  predictions  do  not  appear  to 
correlate  well  with  the  experimental  stress-strain  response  (except  for  very  low  levels  of 
load).  However,  since  only  one  specimen  provided  strain  information,  the  analytical 
predictions  of  stress-strain  response  could  not  be  adequately  compared  with  the  exper¬ 
imental  stress-strain  response. 

Even  though  the  stress-strain  response  plots  (Figures  5-31  through  5-33)  did  not  pro¬ 
vide  as  much  insight  into  the  stress- strain  behavior  of  this  laminate  as  desired,  a  "numer¬ 
ical  instability"  in  the  analysis  is  obvious  in  Figure  5-32.  This  numerical  instability  is 
apparent  as  a  scattering  of  the  data  at  high  levels  of  loading  (Figure  5-32).  These  load 
levels  resulted  in  the  failure  of  the  finite  elements  shown  in  Figure  5-28.  This  numerical 
instability  caused  the  program  to  stop  execution  due  to  the  nonconvergence  of  incremental 
strains  resulting  from  an  addition  increment  of  load  (equation  (41),  Section  n  C).  There¬ 
fore,  for  this  laminate,  the  load  at  which  the  program  stopped  execution  was  considered  to 
be  the  analytical  failure  load. 
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Applied  Stress  (psi) 


Strain  (micro  in/in) 


STRAIN  GAGE  LOCATIONS 
gage#4  (reverse  side) 

W, 

gage#2~ <?>  /• 

gage  #1  — • 

m 

NOTES: 

1.  *  indicates  serial  number  of  test  specimen 

2.  ±45  degree  lay-up  /  eccentricity=0.3"  /  no  mounting  fixture 

3.  no  data  available  from  specimens  034  &  036  (gages  saturated) 


Figure  5-31.  Near  Field  (gage  #1)  Stress  Strain  Response 
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Applied  Stress  (psi) 


Strain  (micro  in/in) 


STRAIN  GAGE  LOCATIONS 
gage#4  (revqse  side) 

m 

Pge#2-a 
gage#l  — •  gage#3 

11 

NOTES: 

1.  *  indicates  serial  number  of  test  specimen 

2.  ±45  degree  lay-up  /  eccentricity=0.3"  /  no  mounting  fixture 

3.  no  data  available  from  specimens  034  &  036  (gages  saturated) 


Figure  5-32.  Near  Field  (gage  #2)  Stress  Strain  Response 
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STRAIN  GAGE  LOCATIONS 
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NOTES: 

1.  *  indicates  serial  number  of  test  specimen 

2.  ±45  degree  lay-up  /  eccentricity=0.3"  /  no  mounting  fixture 

3.  no  data  available  from  specimens  034  &  036  (gages  saturated) 


Figure  5-33.  Far  Field  (gage  #3)  Stress  Strain  Response 
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In  summary,  the  analytical  solution  underestimated  the  failure  loads  (from  16.6%  to 
40.5%)  for  the  [±45°]4S  laminates.  This  underestimation  may  be  due  to  scissoring  of  the 
±45  degree  fibers,  or  may  be  a  result  of  ignoring  the  interlaminar  stresses  (since  this  was  a 
2-D  study).  It  appears  that  the  analytical  solution  predicts  a  shear  (mode  II)  failure  for  this 
stacking  sequence,  and  the  high  strains  at  failure  (necking)  were  also  predicted  analytically. 
Finally,  although  it  appears  that  the  analytical  stress-strain  response  correlated  (for  low 
levels  of  load)  with  the  experimental  stress-strain  response,  the  experimentation  did  not 
provide  enough  useful  data  to  verify  the  analytical  predictions  of  stress-strain  response. 

(4)  [0°/±45o/90°]2s  lanunaics. 

Table  5-10  shows  the  experimental  results  of  the  ultimate  strength  tests  for  the 
[0°/±45o/90°]2s  laminates.  The  scatter  of  the  experimental  data  (Table  5-10)  is  shown 
graphically  in  Figure  5-34.  This  figure  also  lists  the  standard  deviation  for  each  hole 
location  (value  of  eccentricity)  tested  experimentally.  The  standard  deviation  was  calculated 
from  equation  (1),  Section  V  B(l).  Figure  5-34  shows  that  experimental  scatter  was  not 
significant  during  the  testing  of  these  laminates. 

A  summary  of  the  average  experimental  failure  loads  (calculated  from  Table  5-10),  as 
well  as  the  analytical  predictions  of  failure  loads,  are  given  in  Table  5-11.  In  addition. 
Table  5-11  shows  the  percent  difference  between  the  average  experimental  failure  loads  and 
the  analytical  predictions  of  failure  loads.  This  percent  difference  was  calculated  from 
equation  (2),  Section  V  B(l). 

If  the  analytical  predictions  of  the  nominal  (far  field)  failure  stresses  are  examined,  the 
effects  of  eccentricity  on  the  nominal  failure  stress  can  be  determined.  Table  5- 12  lists  both 
the  analytical  failure  loads  and  the  nominal  (far  field)  failure  stresses.  As  previously 
discussed  (Section  V  B(l)),  the  distance  from  the  edge  of  the  discontinuity,  to  the  edge  of 


Table  5-10.  Experimental  Results  for  [0o/±45°/90°]2s  La™*11*165- 

Serial  Number 

Eccentricity 

Mounting 

Fixture 

Failure  Load 
(lbs) 

041 

s 

O 

o 

4250.97 

042 

0.0" 

3944.13 

043 

0.0" 

4008.59 

044 

0.1" 

4832.06 

045 

0.1" 

4846.99 

046 

0.1" 

5445.63 

047 

0.3" 

6683.89 

048 

0.3" 

6617.70 

049 

0.3" 

6750.09 

050 

0.3" 

✓* 

6222.60 

051 

0.3" 

v/ 

6413.33 

052 

0.3" 

v/ 

6719.92 

*  ( \S  indicates  that  the  mounting  fixture  was  used) 

the  specimen  (Figure  4-2),  was  always  greater  than  or  equal  to  the  diameter  of  the  discon¬ 
tinuity  (0.4").  This  distance  was  maintained  so  that  edge  effects  would  not  vary  for  the 
different  values  of  eccentricity  considered  in  this  study.  From  Table  5-12,  note  that  the 
analytical  predictions  of  nominal  failure  stresses  do  not  appreciably  differ  for  the  different 
values  of  eccentricity;  therefore,  since  edge  effects  were  controlled  by  maintaining  at  least 
0.4  inches  between  the  edge  of  the  discontinuity  and  the  edge  of  the  specimen,  it  appears 
that  the  values  of  eccentricity  considered  in  this  study  did  not  appreciably  affect  the  nominal 
failure  stress  for  this  stacking  sequence. 

The  nominal  failure  stress  for  [0°/±45°/90°]2s  laminates,  without  a  discontinuity,  is 
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Figure  5-34.  Experimental  Scatter  for  [0°/±45°  /90 °]  2s Laminates 


Table  5-11.  Summary  of  Average  Experimental  Failure  Loads  and 
Analytical  Failure  Loads  for  [0 °/±45  790°]  ,,  Laminates. 

Eccentricity 

Average  Experimental 
Failure  Load  (lbs) 

Analytical  Failure 
Load  fibs) 

%  diff. 

(Ex.  vs  Anal.) 

0.0" 

(no  fixture) 

4067.90 

3920.27 

+3.6% 

0.1" 

(no  fixture) 

5041.56 

4810.02 

+4.6% 

0.3" 

(no  fixture) 

6683.89 

6388.39 

+4.4% 

0.3" 

(with  fixture) 


6451.95 


6323.25 


+2.0% 


Table  5-12.  Summary  of  Analytical  Failure  Loads  and  Nominal 
(far  field)  Failure  Stresses  for  [0°/±45°/9(f]  ^  Laminates. 

Eccentricity 

Analytical  Failure 

Load  (lbs) 

Nominal  Failure  (far 
field)  Stress  (psi) 

3920.27 

38,891.57 

4810.02 

40,901.53 

0.3" 

(no  fixture) 

6388.39 

42,251.26 

0.3" 

(with  fixture) 

6323.25 

41,820.44 

equal  to  98,400  psi  [17].  Therefore,  for  this  laminate,  the  nominal  failure  stress,  due  to 
the  0.4"  discontinuity,  was  reduced  by  a  factor  of  approximately  2.4  (see  nominal  failure 
stresses.  Table  5-12).  This  is  lower  than  the  reduction  of  nominal  failure  stress  (due  to  a 
circular  discontinuity)  for  an  isotropic  material  of  the  same  geometry;  reported  as  3.46  by 
[15]. 

The  failure  mode  of  these  specimens  was  visible  as  a  shear  (mode  II)  failure  surface 
(Figures  5-35  and  5-36).  However,  the  shear  failure,  mode  n,  for  these  laminates 
(Figures  5-35  and  5-36)  does  not  appear  to  be  as  shear  dominated  as  was  noted  for  the 
[±45°]4S  laminates  (Figures  5-24  and  5-25).  Therefore,  it  appears  that  the  failure  surface 
(region)  is  a  combination  of  the  failure  modes  noted  for  the  [O0^,  [90P16],  and  the 
[±45°]4S  laminates. 

Figures  5-37  through  5-40  show  the  analytical  predictions  of  failure  progression  and 
failure  loads  for  the  [0°ft45o/90°]2s  laminates.  Note,  from  these  figures,  that  the  shear 
(mode  II)  failure  is  not  predicted  analytically.  This  could  be  due  to  several  reasons.  First, 
recall  (Section  V  B(3))  that  the  ±45°  plies  require  large  strains  before  failure;  but,  the  0° 
plies  in  the  [0°/±45°/90°]2s  laminates  prevent  large  strains.  Therefore,  once  the  level  of 
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Figure  5-35.  Failed  [0°/±45°/90°|  ^  Specimens  (no  mounting  fixture) 


Figure  5-36.  Failed  [0°/±45o/90°|  2g  Specimen  (mounting  fixture) 


5-49 


The  shaded  portion  of  the  model  is 
magnified  below,  and  the  displacements 
of  this  se&^on  are  magnified  10  times. 


■  indicates  90  degree  ply  failure  at  3680.20  lbs 
O  indicates  90  degree  ply  failure  at  3920.27  lbs 
C3  indicates  0  degree  ply  failure  at  3920.27  lbs 

NOTE:  Only  the  area  around  the  hole  shown  for  clarity. 
ANALYTICAL  FAILURE  LOAD  WAS  3920.27  LBS. 


Figure  5-37.  Failure  Progression  (quasi  /  e-O.O"  /  no  mounting  fixture) 
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The  shaded  portion  of  the  model  is 
magnified  below,  and  the  displacements 
of  this  section  are  magnified  10  times. 


■  indicates  90  degree  ply  failure  at  4335.02  lbs 
CD  indicates  90  degree  ply  failure  at  4810.02  lbs 
indicates  0  degree  ply  failure  at  4810.02  lbs 

NOTE:  Only  the  area  around  the  hole  shown  for  clarity. 


ANALYTICAL  FAILURE  LOAD  WAS  4810.02  LBS. 


Figure  5-38.  Failure  Progression  (quasi  /  e=0. 1 "  /  no  mounting  fixture) 
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The  shaded  portion  of  the  model  is 
magnified  below,  and  the  displacements 
of  this  section  are  magnified  10  times. 


■  indicates  90  degree  ply  failure  at  6183.87  lbs 
n  indicates  0  degree  ply  failure  at  6183.87  lbs 
®  indicates  90  degree  ply  failure  at  6323.25  lbs 
•  indicates  0  degree  ply  failure  at  6323.25  lbs 

NOTE:  Only  the  area  around  the  hole  shown  for  clarity. 

ANALYTICAL  FAILURE  LOAD  WAS  6323.25  LBS. 


Figure  5-40.  Failure  Progression  (quasi  /  e=0.3"  /  mounting  fixture) 
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loading  is  high  enough  to  cause  failure  of  the  zero  degree  plies  adjacent  (above  and  below) 
to  the  discontinuity,  the  level  of  loading  is  much  higher  than  the  neighboring  ±45°  plies  (or 
90°  plies)  can  sustain.  As  a  result,  Figures  5-37  through  5-40  show  that  the  analytical 
solution  predicted  catastrophic  failure  of  these  laminates  once  the  0°  plies  adjacent  (above 
and  below)  to  the  discontinuity  failed.  Another  possibility  which  might  have  affected  the 
analytical  predictions  of  damage  progression  (Figures  5-37  through  5-40)  is  the  fact  that 
the  two-dimensional  study  did  not  account  for  interlaminar  stresses,  which  for  this  stacking 
sequence,  may  alter  the  failure  mode. 

Stress-strain  response  plots  for  the  eccentricity  =  0.3"  model  were  constructed  (Figures 
5-41  through  5-43).  Figure  5-41  shows  the  experimental  and  analytical  near  field  re¬ 
sponse  at  strain  gage  #1.  The  experimental  and  analytical  near  field  response  at  strain  gage 
#2  is  shown  in  Figure  5-42.  And  finally,  the  experimental  and  analytical  far  field  response 
is  shown  in  Figure  5-43  as  the  average  strains  from  strain  gages  #3  and  #4.  From  these 
three  figures,  note  that  there  is  good,  within  10%,  correlation  (up  to  the  level  of  loading 
which  caused  failure  of  the  0°  plies)  between  the  analytical  predictions  of  axial  stress-strain 
response,  and  the  experimental  axial  stress-strain  response.  The  analytical  predictions  of 
near  field  axial  stress-strain  response  were  found  to  deviate  from  the  experimental  stress- 
strain  response  at  the  level  of  loading  which  corresponded  to  the  analytical  failure  of  the 
zero  degree  plies  adjacent,  above  and  below,  to  the  discontinuity  (0°  ply  failures  are  shown 
in  Figure  5-39).  This  deviation  is  shown  in  Figures  5-41  and  5-42.  Finally,  there  was  ex¬ 
cellent  correlation  between  the  analytical  predictions  of  transverse  stress-strain  response 
and  the  experimental  stress-strain  response. 

In  summary,  the  analytical  solution  did  a  good  job  of  predicting  the  failure  loads  for  the 
(0°/t45°/90P]2*  laminates.  The  analytical  stress-strain  response  predictions  correlated  well 
with  the  experimental  stress-strain  response  up  to  the  point  of  0°  ply  failures,  adjacent  to 
the  discontinuity. 
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Figure  5-41.  Near  Field  (gage  #1)  Stress  Strain  Response 
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|  Figure  5-42.  Near  Field  (gage  #2)  Stress  Strain  Response 
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Figure  5-43.  Far  Field  (gage  #3)  Stress  Strain  Response 
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(C)  Eccentricitv/Mounting  Fixture  -  -  Effects. 
(1)  Eccentricity  Effects. 


As  discussed  in  detail  in  Section  m  E  (1),  the  eccentric  circular  discontinuity  creates  a 
centroid  (area)  shift.  This  shift  in  centroid  is  shown  in  Figure  3-11. 

The  eccentricity  effects  and  resulting  bending  moment  (Figure  3-11),  created  higher 
tensile  stresses  above  the  discontinuity.  Therefore,  elements  above  the  discontinuity  tended 
to  fail  more  quickly  than  corresponding  elements  below  the  discontinuity.  These  effects 
are  quite  obvious  from  the  failure  progression  plots  (Figures  5-16, 5-17, 5-18, 5-27,  5-28, 
5-29,  5-39,  5-40). 

From  the  experimentation,  the  eccentricity  effects  (in  terms  of  failure  progression)  were 
not  observed.  This  was  due  to  the  fact  that  failure  was  too  rapid  to  determine  if  (by 
viewing  the  video  tape  frame  by  frame)  failure  initiated  above  or  below  the  discontinuity. 
Only  a  limited  quantity  of  Gr/PEEK  panels  were  available  for  testing.  Therefore,  all 
specimens  were  tested  to  their  ultimate  failure  load.  If  more  Gr/PEEK  panels  had  been 
available  during  this  study,  the  experimental  verification  of  the  analytical  predictions  of 
eccentricity  effects  on  failure  progression  may  have  been  possible.  This  could  have  been 
achieved  by  testing  specimens  to  a  percentage  of  their  ultimate  failure  load.  These 
specimens  could  have  then  been  examined  for  damaged  areas  through  the  use  of  stereo  x- 
ray  techniques  [10,13],  and  failure  initiation,  as  well  as  failure  progression,  could  have 
been  determined. 

(2)  Mounting  Fixture  Effects. 

The  linear  study.  Section  III  E,  predicted  that  the  use  of  the  mounting  fixture  would  not 
appreciably  affect  the  stress  state  near  the  discontinuity.  However,  the  nonlinear  analysis 
predicted  slighdy  lower  failure  loads  for  models  which  utilized  the  boundary  conditions 
(Section  HI  D)  associated  with  the  mounting  fixture.  These  analytical  failure  loads  are 
summarized  (from  Tables  5-2, 5-5, 5-8,  and  5-1 1)  in  Table  5-13.  These  lower  failure  load 
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Table  5-13.  Summary  of  Analytical  Predictions  of  Failure  Loads 
with  and  without  die  Mounting  fixture. 

Stacking 

Sequence 

Analytical  Failure  Load 
(lbs)  for  eccentricity  of 
0.3"  (mounting  fixture) 

Analytical  Failure  Load 
(lbs)  for  eccentricity  of 

0.3"  (no  mounting  fixture) 

8952.14 

9241.27 

[90°16  3 

1172.56 

1217.37 

[±45°  ] 

4s 

2590.33 

2931.39 

[0°/±45o/90°] 

2$ 

6323.25 

6388.39 

predictions  (for  models  which  used  the  mounting  fixture  boundary  conditions)  can  be 
explained  in  terms  of  the  in-plane  rotation  which  was  permitted  by  the  mounting  fixture. 
The  output  data  (for  all  four  stacking  sequences)  was  examined,  and  depending  on  the 
boundary  conditions  used  (i.e.  in-plane  rotation  was/was  not  permitted)  a  difference  in  near 
field  axial  stresses  was  noted.  These  differences  in  axial  stresses  can  be  explained  by 
considering  Figure  5-44.  This  figure  shows  two  [±45°]4S  models  (at  the  same  level  of 
load).  One  of  these  models  (Figure  5-44  A)  uses  the  boundary  conditions  associated  with 
the  mounting  fixture,  while  the  other  model  (Figure  5-44  B)  uses  the  boundary  conditions 
associated  with  no  mounting  fixture.  In  Figure  5-44  A,  note  that  in-plane  rotation  is  very 
apparent,  while  Figure  5-44  B  shows  that  with  the  no  mounting  fixture  boundary 
conditions,  no  in-plane  rotation  is  allowed  In  addition,  after  examining  the  stress  levels, 
for  both  models,  above  and  below  the  discontinuity  (Figure  5-44  C),  the  axial  stresses  for 
the  model  which  allowed  in-plane  rotation  were  found  to  be  higher,  above  the  discon- 
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Figure  5-44  A.  [±45° 1  Model  (eccentricity^O  3" /mounting  fixture) 

4s 

initial  displ  =  .002866"/displacements  magnified  1000  times/NOTE  ROTATION 


Figure  5-44  B.  [±45°]  Model  (eccentricity=03"lno  mounting  fixture) 

4s 

initial  displ  =  .002866'Vdisplacements  magnified  1000  times/NO  ROTATION 


Figure  5-44  C.  Elements  Adjacent  to  Discontinuity  (undistorted) 


Figure  5-44.  Mounting  Fixture  Effects 
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tinuity,  than  for  the  model  which  did  not  allow  in-plane  rotation .  Furthermore,  the  axial 
stresses  for  the  model  which  allowed  in-plane  rotation  were  found  to  be  lower,  below  the 
discontinuity,  than  for  the  model  which  did  not  allow  in-plane  rotation.  Therefore,  the  in¬ 
plane  rotation  creates  higher  axial  stresses  above  the  discontinuity  and  consequently,  leads 
to  failure  at  lower  levels  of  load  than  for  the  models  which  do  not  allow  in-plane  rotation 
(see  Table  5-13).  As  noted  in  Figure  5-44,  the  displacements  for  Figures  5-44  A  and  5-44 
B  were  magnified  1000  times.  If  the  displacements  were  not  significantly  magnified  (over 
100  times),  the  in-plane  rotation  was  not  visible. 

Experimentally,  the  use  of  the  mounting  fixture  did  not  always  lower  the  failure  loads 
for  the  laminates  (see  Tables  5-2, 5-5, 5-8,  and  5-11).  This  is  most  likey  due  to  the  fact 
that  only  three  specimens  were  tested  experimentally  fra-  each  of  the  stacking  sequences 
(see  Tables  5-1, 5-4, 5-7, 5-10).  Therefore,  since  the  analytical  solution  predicted  only 
slightly  lower  failure  loads  (with  the  use  of  the  mounting  fixture),  even  small  scattering  of 
the  experimental  data  (see  Figures  5-5, 5-12, 5-22,  and  5-34)  could  account  for  not  always 
noting  the  trend  of  lower  failure  loads  with  the  use  of  the  mounting  fixture.  Also,  as  noted 
in  the  previous  paragraph,  the  in-plane  rotations  were  very  small,  and  therefore,  it  was 
impossible  to  observe  if  in-plane  rotations  actually  occurred  during  the  experimentation. 
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VI.  Conclusions 


This  section  summarizes  the  conclusions  relating  to  the  objectives  listed  in  Section  I  A. 
Hrst,  the  conclusions  from  a  comparison  of  analytical  and  experimental  results  are  pre¬ 
sented.  Secondly,  conclusions  relating  to  both  boundary  condition  effects  (Le.  mounting 
fixture  effects)  and  eccentricity  effects  on  the  failure  characteristics  of  the  Gr/PEEK  lami¬ 
nates  are  discussed.  And  finally,  general  conclusions  relating  to  die  suitability  of  the 
numerical  tools  (linear  and  nonlinear  finite  element  programs),  used  in  this  study,  are 
discussed. 

A.  Experimental  versus  Analytical  Failure  Predictions. 

(1)  rO°i6l  laminates. 

The  finite  element  (FE)  solution  did  a  good  job  of  predicting  the  splitting  (failure)  loads 
for  these  laminates.  In  addition,  the  FE  stress-strain  response  predictions  correlated  well 
with  the  experimental  stress-strain  response.  However,  the  FE  solution  did  not  predict  the 
failure  progression  (Le.  splitting)  in  terms  of  damage  accumulation.  The  splitting  of  these 
laminates  may  be  caused  by  a  dynamic  phenomenon,  such  as  a  stress  wave,  or  may  be 
caused  due  to  some  other  three-dimensional  effect,  such  as  edge  effects.  Therefore,  the 
two-dimensional  FE  solution  was  not  able  to  predict  failure  progression. 

(2)  [90°  j  5]  laminates. 

For  these  laminates,  the  FE  solution  did  a  good  job  of  predicting  both  the  experimental 
failure  loads  and  the  experimental  stress-strain  response.  In  addition,  the  numerical  pre¬ 
dictions  of  damage  accumulation  (failure  progression)  showed  that  damage  initiation  would 
occur  above  and  below  the  discontinuity,  and  that  failure  would  occur  rapidly  once  die 
damage  initiation  occurred.  These  predictions  (of  failure  progression)  matched  the  results 
obtained  from  experimentation,  where  is  was  observed  that  these  laminates  failed  suddenly, 
and  damage  initiation  occurred  above  and  below  die  discontinuity. 
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(3)  f±45%  laminates. 

For  these  laminates,  die  FE  solution  was  found  to  significantly  underestimate  the 
experimental  failure  loads.  This  underestimation  may  be  due  to  scissoring  of  the  ±45 
degree  fibers,  or  may  be  a  result  of  ignoring  the  interlaminar  stresses  (since  this  was  a  two- 
dimensional  study).  Since  the  scissoring  phenomenon  causes  die  fibers  to  become  more 
aligned  with  the  axial  loading,  one  would  expect  that  die  experimental  failure  loads  would 
be  higher  than  the  analytical  predictions  (since  the  program  utilized  did  not  account  for 
scissoring).  The  experimental  scattering  of  test  data  was  not  significant  enough  to  account 
for  the  numerical  underestimation  of  the  failure  loads. 

The  FE  solution  predicted  a  shear  (mode  II)  failure  for  these  laminates,  and  also  pre¬ 
dicted  the  high  strains  (necking)  at  failure.  Finally,  although  the  FE  predictions  of  stress- 
strain  response  correlated,  for  low  levels  of  load,  with  the  experimental  stress-strain 
response,  the  experimental  on  did  not  provide  enough  useful  data  to  verify  die  analytical 
predictions  of  stress-strain  response. 

(4)  [0°/t45°/90^2s  laminates. 

The  FE  solution  did  a  good  job  of  predicting  the  experimental  failure  loads  for  these 
laminates.  In  addition,  the  FE  stress-strain  response  predictions,  up  to  the  point  of  failure, 
correlated  well  with  the  experimental  stress-strain  response. 

The  failure  of  these  laminates  was  characterized  by  a  shear  (mode  II)  failure.  However, 
die  failure  regions  of  these  laminates  did  not  appear  to  be  as  shear  dominated  as  was  noted 
for  die  (±45°]4I  humnates.  Therefore,  it  appears  that  the  failure  mode  for  these  laminates  is 
some  combination  of  the  failure  modes  noted  for  die  [0°16],  [90°16],  and  [±45°]4S 
laminates. 

The  FE  solution  did  not  predict  the  shear  (mode  II)  failure.  This  may  be  due  to  several 
reasons.  First,  the  ±45  degree  plies  require  large  strains  before  failure;  but,  the  zero  degree 


plies  in  the  [0P#45°/90P)2s  laminate  prevent  large  strains.  Therefore,  once  the  level  of 
loading  was  high  enough  to  cause  failure  of  die  zoo  degree  plies,  adjacent  (above  and 
below)  to  the  discontinuity,  the  level  of  load  was  much  higher  than  the  neighboring  ±45 
degree  plies  (or  90  degree  plies)  could  sustain.  As  a  result,  the  FE  solution  predicted  cat¬ 
astrophic  failure  of  these  laminates  once  the  zero  degree  plies,  adjacent  (above  and  below) 
to  the  discontinuity,  failed.  Another  possibility  which  might  have  affected  the  FE  pre¬ 
dictions  of  damage  progression,  is  the  fact  that  the  two-dimensional  FE  solution  did  not 
account  for  interlaminar  stresses,  which  for  this  stacking  sequence,  may  alter  the  failure 
mode. 

During  the  experimentation,  audible  "popping"  of  die  fibers  could  be  heard  for  several 
minutes  before  the  laminates  failed.  Therefore,  one  would  think  that  the  damage  zone  (at 
die  failure  load)  would  not  be  as  small  as  the  damage  zone  predicted  by  die  FE  solution. 

One  possibility  for  the  prediction  of  a  small  damage  zone,  is  the  fact  that  interlaminar 
stresses  were  ignored  in  this  study,  and  as  mentioned  previously,  these  stresses  may  have 
an  effect  on  the  failure  characteristics  of  this  stacking  sequence. 

b.  Boundary  Condition/EcangicitY  Efforts- 

(1)  Boundary  Condition  Effects. 

The  linear  FE  study  predicted  that  the  use  of  the  mounting  fixture  would  not  appreciably 
affect  the  stress  states  (near  the  discontinuity)  for  the  Gr/PEEK  laminates  considered  in  this 
study.  The  nonlinear  FE  study  predicted  that  slightly  lower  failure  loads  would  result 
from  die  use  of  die  mounting  fixture.  These  lower  failure  loads  appear  to  result  from  the 
in-plane  rotation  allowed  by  die  mounting  fixture.  This  in-plane  rotation  results  in  a  higher 
axial  stress  region  above  die  discontinuity,  and  consequently  led  to  failure  at  lower  levels  of 
load. 
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The  experimental  verification  of  the  FE  predictions  of  boundary  condition  effects,  was 
inconclusive.  Since  only  three  specimens  (for  each  stacking  sequence  and  value  of  ecc¬ 
entricity)  were  tested  experimentally,  experimental  scattering  of  test  data  could  account  for 
not  always  observing  a  lowering  of  the  experimental  failure  loads  through  die  use  of  the 
mounting  fixture.  It  should  be  noted  that  the  experimentation  did  verify  lower  failure  loads 
for  die  [0°|$]  and  (Q^>/k4So/90°]2s  laminates  when  die  mounting  fixture  was  used. 

(2)  Eccentricity  Effects. 

The  eccentric  circular  discontinuity  created  a  centroid  (area)  shift  This  centroid  shift 
and  resulting  bending  moment,  created  higher  tensile  stresses  above  die  discontinuity.  This 
phenomenon  was  predicted  by  the  FE  study,  and  it  was  observed  that  elements  above  the 
discontinuity  tended  to  fail  mare  quickly  than  corresponding  elements  below  the  discon¬ 
tinuity.  Therefore,  eccentricity  effects  led  to  a  higher  axial  stress  region  above  the  dis¬ 
continuity.  However,  the  nominal  (far  field)  failure  stresses,  for  all  four  stacking 
sequences,  did  not  apprc.  *  ’y  differ  for  the  different  values  of  eccentricity  considered. 

Eccentricity  effects  (in  terms  of  failure  progression)  were  not  observed  experimentally. 
This  was  due  to  the  fact  that  failure  was  too  rapid  to  determine  if  (by  viewing  a  video  tape 
of  the  testing  frame-by-frame)  failure  initiated  above  or  below  the  discontinuity.  If  more 
panels  of  Gr/PEEK  had  been  available  for  this  study,  stereo  x-ray  techniques  [10,13]  could 
have  been  used  to  study  specimens  which  had  been  tested  to  a  percentage  of  their  ultimate 
load.  This  may  have  provided  insight  into  die  eccentricity  effects  on  the  initiation  and 
progression  of  failure. 

c.  General  Conclusions- 
(1)  linear  Haitt  Element  Promm  (PLSTR). 

The  linear  finite  element  program  was  very  useful  in  providing  a  numerical  tool  which 
was  used  to  develop  the  finite  element  models  and  predict  states  of  stress  for  this  study.  In 
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addition,  the  program  was  used  to  generate  data  for  contour  plotting  and  deformed  model 
plots  (models  under  loading).  Since  inexpensive  computer  resources  (such  as  VAX  main¬ 
frames)  can  quickly,  under  60  real  time  seconds,  handle  the  computations  required  by  this 
program,  the  program  is  very  useful  for  preliminary  work  leading  into  a  nonlinear  study 
(where  more  expensive  computer  resources  such  as  the  CRAY  are  required).  Finally,  the 
linear  finite  element  program  closely  approximated  the  experimentally  obtained  stress  states 
of  the  Gr/PEEK  laminates. 

(2)  Nonlinear  Material  Finite  Element  Program  (PLSTREN). 

The  nonlinear  material  finite  element  program  is  a  progressive-ply-failure  program 
which  predicts  both  damage  initiation  and  progression.  This  program  was  found  to  give 
good  results  for  predicting  the  failure  characteristics  (due  to  an  eccentric  circular  discon¬ 
tinuity)  for  [0P16],  [90^ 16],  and  [0°/±45°/90P]2s  Gr/PEEK  laminates.  Since  this  program 
requires  efficient  computer  resources,  a  detailed  linear  study  (using  PLSTR)  should  be 
conducted  before  using  PLSTREN.  By  accomplishing  this  preliminary  linear  work,  the 
lengthy  modeling  process,  and  inevitable  errors,  can  be  completed  on  low  cost  computer 
resources.  Finally,  the  total  strain  energy  failure  criterion  and  the  unloading  option  utilized 
by  PLSTREN  produced  relatively  good  results. 

(3)  Closing. 

This  thesis  has  taken  a  detailed  "look"  at  the  objectives  stated  in  Section  I  A.  It  has 
shown  that  the  nonlinear  material  finite  element  program  (PLSTREN)  is  able  to  predict  the 
eccentricity  effects  on  the  failure  characteristics  of  Gr/PEEK  laminates.  It  is  hoped  that  this 
study  will  contribute,  in  some  small  way,  to  other  research  projects  involving  thermoplastic 
composite  laminates,  such  as  Gr/PEEK. 
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Appendix  A.  Material  Property  Curves 


Basic  property  tests  were  conducted  by  [13]  to  obtain  the  material  properties  for 
Gr/PEEK  (at  room  temperature).  This  data  is  shown  graphically  (Figures  A1  through  A7) 
to  demonstrate  the  nonlinear  relationships  between  stress  and  strain  for  Gr/PEEK. 


Figure  A-l.  Oj  vs  (0°  tension) 


A-l 


Transverse  Strain  (in/in) 
Figure  A-4.  a2  vs  (90 0  compression) 


Figure  A-5.  x12  vs  y12  (±45  tension) 


Poisson's  Ratio 


Figure  A-6.  v12  vs  Ej  (0°  tension) 


Strain  (in/in) 

Figure  A-7.  v12  vs  (0°  compression) 
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Appendix  B.  Contour  Plots  (Linear  Analysis) 


As  mentioned  in  Section  in  E  (1),  contour  plots  were  generated  for  all  of  the  specimens 
tested  experimentally  (see  Tables  4-1  through  4-4).  These  contour  plots  were  generated  to 
study  the  effects  of  both  the  eccentric  circular  discontinuity  and  the  mounting  fixture  on  the 
overall  stress  state  of  the  Gr/PEEK  laminates.  Since  the  highest  concentration  of  stress 
was  located  adjacent  to  the  discontinuity,  a  window  (see  Figure  B-l)  surrounding  the 
discontinuity  was  selected  for  the  contour  plots.  The  contour  plots  (see  Figures  B-2 
through  B- 1 1 )  are  grouped  according  to  the  stacking  sequence  of  the  laminate  and  also  by 
the  boundary  conditions  used  with  the  finite  element  models.  These  boundary  conditions 
simulated  the  use  of  either  the  mounting  fixture  or  the  Instron  Test  Machine  grips  (see 
Section  m  D).  The  term  "no  hinge"  on  the  contour  plots  shown  in  this  section  indicates 
that  no  in-plane  rotation  (analytically)  was  permitted.  The  term  "hinge"  on  the  contour 
plots  indicates  that  ;n-plane  rotation  (analytically)  was  permitted.  Finally,  the  value  of 
eccentricity  is  also  given  on  each  contour  plot 


Window  for  contour  plots 


Window  for  contour  plots 


Window  for  contour  plots 


Eccentricity  =  0.3"  Model 


Figure  B-5.  Contour  Plots  (90  degree  lay-up /mounting  fixture) 


Figure  B-6.  Contour  Plots  (±45  degree  lay-up  /  no  mounting  fixture) 
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Figure  B-8.  Contour  Plots  (quasi  lay-up  /no  mounting  fixture) 


Figure  B-9.  Contour  Plots  (quasi  lay-up  /  mounting  fixture) 
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Appendix  C.  Stress  State  Plots  (Linear  Analysis) 


As  mentioned  in  Section  m  E  (2)a,  plots  were  generated  for  both  near  field  and  far  field 
stress  states.  These  plots  were  made  in  order  to  study  boundary  condition  and  eccentri¬ 
city  effects  near  the  discontinuity  (near  field)  as  well  as  near  die  tab  area  (far  field).  Figure 
C-l  shows  the  near  field  location.  Figures  C-3  through  C-6  were  plotted  at  the  near  field 
location  shown  in  Figure  C-l.  Figure  C-2  shows  the  far  field  location.  Figures  C-7 
through  C-10  were  plotted  at  the  far  field  location  shown  in  Figure  C-2.  The  data  used  to 
generate  Figures  C-3  through  C-10,  was  obtained  from  output  data  from  linear  finite 
element  runs. 
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Figures  C-3  through  C-6  were  plotted  from  data  obtained  along 
the  line  X=7.0”. 


Figure  C- 1 .  Near  Field  Locations  for  Plots 


Figure  C-2.  Far  Field  Locations  for  Plots 
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Figure  C-3.  Near  Field  Axial  Strain  Plots  (0  degree  lay-up) 
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Figure  C-8.  Far  Field  Axial  Strain  Plots  (90  degree  lay-up) 
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Appendix  D.  Output  Files  from  Nonlinear  Runs 


This  appendix  discusses  the  output  files  created  during  a  nonlinear  finite  element  run. 
Once  a  nonlinear  run  had  been  completed  by  the  CRAY,  the  following  output  files  were 
created. 

(1)  Tape  6  (output  data) 

(2)  Tape  3  (used  to  restart  a  run  from  the  point  of  termination) 

(3)  Tape  9  (nodal  displacement  data  as  a  function  of  the  load  increment) 

(4)  Listing  (source  code  listing  and  any  error  messages) 

(1)  Tape  6. 

The  output  file  (Tape  6)  contained  the  following  information. 

a) .  Summary  of  input  data 

aa).  control  card  summary 
ob).  x  &  y  location  of  nodes 
cc).  boundary  conditions 

dd).  connectivity  (nodes  surrounding  the  element) 
ee).  material  properties 

b) .  For  each  load  increment  a  summary  of 

aa).  displacement  of  each  node 

bb).  load  tor  each  node 

cc).  global  and  local  stresses  (each  element) 

dd).  global  and  local  strains  (each  element) 

ee).  strain  energy  (each  element) 

ft),  element  failure  message  (if  applicable) 

(2)  Tape  3. 

An  execution  time  limit  (prescribed  by  the  user)  was  incorporated  in  the  batch  file.  This 
time  limit  was  used  to  prevent  long  execution  times  for  input  files  containing  errors.  If  this 
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time  limit  was  exceeded,  execution  would  be  stopped  and  a  Tape  3  produced.  This  tape 
could  then  be  used  to  restart  the  program  from  the  point  of  termination.  This  "restart" 
capability  saved  money  (by  saving  computer  time)  for  runs  in  which  the  modeled  laminate 
had  not  failed  at  the  time  the  program  was  stepped. 

(3)  Tape  9. 

Tape  9  contained  the  nodal  displacements  of  each  node  in  the  finite  element  model  (for 
each  increment  of  loading).  This  file  was  used  for  making  distorted  plots  of  the  model  at 
any  increment  of  loading.  After  a  nonlinear  run,  this  file  was  transferred  to  a  VAX  main¬ 
frame  computer,  and  was  used  as  the  input  data  file  to  a  plotting  routine  (written  by  Dr. 
Sandhu).  Section  V  B  contains  distorted  plots  (for  the  various  stacking  sequences)  which 
were  produced  from  Tape  9  data. 

(4)  Listing. 

The  last  output  file  was  the  listing.  This  file  contained  a  listing  of  the  source  code  for 
the  nonlinear  program.  In  addition,  this  file  contained  the  sequence  of  steps  executed  by  _ 
the  CRAY  during  the  run  process.  This  listing  was  an  invaluable  source  of  information 
(especially  to  this  student)  for  "debugging"  purposes. 
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Appendix  E.  Theory  of  Photoelasticitv 


Stress-Optics  (ohotoclasticitv). 

Photoelasticity  can  be  considered  an  experimental  stress  analysis  technique.  The 
versatility  of  photoelastic  coatings  was  utilized  in  this  thesis  to  provide  a  visible  picture  of 
surface  stress  distribution.  Photoelastic  fringe  patterns  were  compared  to  contour  plots 
produced  from  the  analytical  solution  (see  Section  V  A). 

Many  materials,  notably  plastics,  are  optically  isotropic  when  they  are  unstressed  but 
become  optically  anisotropic  when  stressed.  This  results  in  a  change  in  the  index  of 
refraction.  Therefore,  the  index  of  refraction  becomes  a  function  of  the  applied  load.  In 
short,  this  becomes  die  underlying  principal  of  the  theory  of  photoelasticity.  When  a 
polarized  beam  of  light  passes  through  a  photoelastic  coating,  the  beam  splits  and  two 
polarized  beams  are  propagated  in  the  direction  of  the  principal  strains. 

Brewster's  Law  states  that  "The  relative  change  of  refraction  is  proportional  to  the 
difference  of  principal  strains."  Based  upon  die  strain  intensities  (in  the  direction  of 
principal  strains),  and  the  speed  of  light  vibrating  in  these  directions,  the  time  necessary  for 
each  of  the  polarized  beams  to  pass  through  the  photoelastic  coating  will  be  different  The 
relative  retardation  between  these  beams  causes  interference  and  results  in  visible  fringes. 
Therefore,  although  photoelastic  coatings  provide  a  visible  picture  of  surface  stress  distri¬ 
bution,  the  actual  fringes  observed  during  photoelastic  testing  are  created  by  the  difference 
in  principal  strains.  This  retardation  can  be  calculated  as: 

5 = c  [v7 '  V7] = t(n*  *  ny} 

where: 

8  ■  retardation 

C«  speed  of  light 

t »  photoelastic  coating  thickness 

Vx’  ■  speed  of  light  in  direction  of  the  first 
principal  strain 
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Vy'  =  speed  of  light  in  direction  of  the  second 
principal  strain 
C 

nx  =  V7 

C 

ny=V7 

Furthermore,  Brewster’s  Law  can  be  stated  mathematically  as: 

(nx-ny)  =  k(Ex'-ey') 

where: 

k  =  strain-optical  coefficient 

Therefore,  the  expression  commonly  used  for  strain  measurements  is  given  by: 

ex'-ey’=^E 

where: 

Ex'  *  first  principal  strain 
Ey'  ~  second  principal  strain 
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Since  eccentricity  effects  on  the  failure  modes  and  failure  loads  of  the  Gr/PEEK 
laminates  were  of  interest  in  this  study,  the  analytical  study  attempted  to  predict  the 
behavior  of  the  Gr/PEEK  laminates  due  to  an  eccentric  circular  discontinuity.  The  eccentric 
circular  discontinuity  created  a  centroid  (area)  shift  in  the  Gr/PEEK  laminates.  This  shift  in 
centroid,  combined  with  the  axial  loading,  created  a  bending  moment  which  varied 
depending  on  the  eccentricity  of  the  circular  discontinuity.  The  analytical  study  was  able  to 
accurately  predict  the  failure  modes  and  failure  loads  resulting  from  eccentricity  effects. 
Although  the  final  failure  mode,  and  failure  load,  was  available  from  experimentation,  the 
failure  progression  resulting  from  eccentricity  effects  was  not  observed  during  exper¬ 
imentation.  All  experimentation  was  recorded  on  video  tape  for  post  failure  analysis. 
However,  failure  was  too  rapid  to  observe  failure  progression  (by  viewing  the  video  tape 
frame-by-frame). 


